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INTRODUCTION 


The  formalities  associated  with  this  contract  were  completed  during  October  1986, 
and  the  effective  starting  date  of  this  project  is  October  1,  1986.  Recently,  we  requested, 
and  were  granted  a  no-cost  extension  until  December  31,  1989.  The  present  final  report 
covers  the  period  from  the  starting  date  until  October  1,  1986  to  December  15,  1989. 

The  objectives  of  this  project  were: 

•  to  ascertain  whether  the  course  of  spontaneous  or  phytohemagglutinin  (PHA)-induced 
blastic  transformation  of  human  lymphocytes  in  vitro  is  affected  by  exposure  to  con¬ 
tinuous  wave  (CW)  or  pulsed  (PW)  2450  MHz  microwaves; 

•  to  compare  effects  obtained  by  CW  and  PW  microwave  exposure  at  the  same  average 
specific  absorption  rate  (SAR); 

•  to  compare  effects  obtained  by  microwave  heating  and  conventional  heating  needed 
to  increase  the  temperature  of  the  sample  by  0.5,  1.0,  1.5  and  2.0°C. 

•  to  ascertain  whether  chromosomal  aberrations  are  induced  under  the  conditions  of 
these  experiments. 

In  order  to  provide  reliable  quantitative  data  on  exposure,  a  system  with  automated 
dosimetry  was  developed,  and  tested  for  biocompatibility  with  human  lymphocyte  cultures. 
A  method  for  quantitation  of  lymphocyte  transformation  based  on  image  analysis  was 
developed  and  tested.  Experiments  on  effects  of  CW  or  PW  2450  MHz  exposure  at  heating 
and  non-heating  levels  were  carried  out.  Effects  of  conventional  heating  were  examined 
and  compared  with  effects  of  microwave  heating.  For  convenience  and  clarity  the  work 
performed  is  described  in  the  following  sequence:  (1)  exposure  system  and  dosimetry,  and 
(2)  studies  on  lymphocyte  cultures,  and  (3)  conclusions. 

EXPOSURE  SYSTEM  AND  DOSIMETRY 
Description  of  the  Exposure  System 

The  experiments  planned  in  this  project  necessitated  the  design  and  assembly  of  an 
exposure  system  that  would  meet  several  engineering,  dosimetric  and  biological  require¬ 
ments.  The  system  should  provide: 

(1)  an  environment  compatible  with  lymphocyte  culture  in  vitro  in  respect  to  relative 
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humidity  (90-100%),  CO2  content  in  air  (5%),  and  temperature  (37°C),  as  well  as  the 
possibility  of  heating  and  equilibrating  the  temperature  of  the  samples  up  to  39° C 
using  microwave  or  conventional  heating; 

(2)  the  possibility  of  exposing  samples  (preferably  involving  the  exposure  of  multiple  sam¬ 
ples  under  identical  conditions)  to  CW,  or  PW  2450  MHz  with  a  high  ratio  of  peak 
to  average  power,  and  the  possibility  of  simultaneous  sham  exposure  of  control  sam¬ 
ples  under  strictly  comparable  conditions,  where  the  only  difference  between  exposed 
and  sham  exposed  samples  would  consist  in  the  presence  (or  absence)  of  microwave 
absorption  in  the  sample; 

(3)  continuous  measurements  and  registration  of  the  temperature  of  exposed  and  sham 
exposed  samples  during  exposure  to  provide  a  temperature  profile  over  time  history 
of  the  samples; 

(4)  continuous  monitoring  of  the  wave  form  and  repetition  rate  in  the  case  of  PW  expo¬ 
sures. 

The  first  and  second  requirements  were  met  by  placing  two  shorted  rectangular  S-band 
waveguides  300  mm  long  (Omega  Laboratories  Models  6101)  in  a  conventional  CO2  tissue 
culture  incubator  (Forma  Scientific)  thermostated  at  37° C.  A  styrofoam  block  placed  in 
the  waveguide  serves  as  support  for  the  sample  holder,  and  ensures  that  the  sample  is 
centrally  placed  in  the  waveguide  in  the  same  location  in  successive  experiments.  A  small 
hole  drilled  in  the  waveguide  allows  the  introduction  of  a  minimally-perturbing  temperature 
probe.  The  hole  is  drilled  near  the  guide  wall,  at  the  center  of  the  shorter  dimension,  this  is 
a  low  electric  field  point  and  the  hole  produces  minimal  disturbance.  One  of  the  waveguides 
serves  for  sham  exposures,  the  other  one  is  connected  through  a  matched  coaxial  cable  to 
the  microwave  power  source.  This  consists  of  several  elements.  The  output  from  a  CW 
Hewlett-Packard  Model  8616A  oscillator  feeds  a  Hughes  Model  1177H  travelling  wave  tube 
amplifier.  Pulsed  or  amplitude  modulated  waveforms  can  be  obtained  through  the  use  of 
a  Hewlett-Packard  Model  8403A  p-i-n  modulator,  for  high  peak  power  pulsed  exposures, 
the  signal  can  be  further  amplified  with  a  1  kW  amplifier  (MCL,  La  Grange,  IL).  The 
waveguide  receives  its  power  through  an  isolator  and  a  dual  directional  coupler  that  allows 
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forward  and  reflected  powers  to  be  measured  with  Hewlett-Packard  Model  432A  power 
meters  and  waveforms  to  be  monitored  using  a  Phillips  PM  3240  oscilloscope.  The  amount 
of  power  supplied  to  the  waveguide  is  controlled  by  the  output  from  the  oscillator  and  the 
voltage  of  the  1  kW  amplifier,  and  monitored  by  measurements  of  the  forward  and  reflected 
powers,  and  by  the  readings  on  the  oscilloscope.  The  main  function  of  the  oscilloscope  is, 
however,  to  monitor  waveforms.  These  means  of  control  and  monitoring  are  only  auxiliary, 
as  the  conditions  of  exposure  are  characterized  by  dosimetry  based  on  measurements  of 
temperature  and  of  SAR  within  the  sample,  as  specified  by  the  3rd  requirement  the  system 
should  meet. 

The  dosimetric  part  of  the  system  consists  of  two  minimally-perturbing  temperature 
probes,  connected  to  a  Hewlett-Packard  Model  59306A  relay  activator,  which  is  itself  con¬ 
nected  to  a  Keithley  Model  192  digital  voltmeter  (DVM).  Both  the  relay  activator  and  the 
DVM  are  under  the  control  of  a  Hewlett-Packard  Model  86  desk-top  computer  through  the 
HPIB  (IEEE-488)  bus.  Initially  a  Narda  Model  8011B  non-perturbing  double  temperature 
probe  was  used,  later  two  Vitek  Model  101  probes  were  substituted  for  the  Narda  probe. 
The  Narda  probe  requires  cumbersome  calibration  against  a  NBS  standard  and  introduc¬ 
tion  of  correction  factors  to  relate  the  current  to  temperature,  a  linear  relationship  being 
maintained  over  limited  temperature  ranges.  Vitek  probes  are  much  easier  to  calibrate, 
are  more  stable,  and  the  current-temperature  relationship  is  linear  over  the  temperature 
range  of  interest. 

Under  computer  control  the  temperatures  of  exposed  and  sham  exposed  samples  are 
recorded  sequentially  and  stored  in  memory  (records  are  maintained  by  keeping  a  library 
of  micro-floppy  disks).  Since  the  temperature  of  the  sham  exposed  sample  is  essentially 
constant  in  view  of  the  stability  of  the  temperature  of  the  conventional  incubator,  it  can  be 
monitored  less  frequently  than  the  temperature  of  the  exposed  samples.  For  conventional 
thermal  exposures  the  temperature  can  be  elevated  by  increasing  the  temperature  of  the 
incubator,  maintaining  control  samples  in  another  CO]  incubator. 

The  sample  holder  has  to  be  a  tissue  culture  vessel  made  of  materials  non-toxic  to 
cells  in  vitro  and  transparent  to  2450  MHz  microwaves.  The  geometry  of  the  container  has 
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to  allow  for  uniform  energy  deposition  within  the  sample.  These  considerations  limited 
the  choice  of  possible  sample  holders.  Previous  experience  with  tissue  culture  exposures  in 
a  S-band  waveguide  demonstrated  that  the  use  of  T-flasks  or  Petri  dishes  leads  to  signifi¬ 
cant  nonuniformities  in  energy  deposition.  Two  types  of  sample  holders  were  tested,  both 
rectangular  tissue  culture  Lab-Tek  (Miles  Laboratories)  chamber  slides,  a  two  chamber 
model  No4802,  and  a  four  chamber  model  No4804.  Analysis  of  SAR  distribution  based 
on  measurements  of  temperature  (see  below)  in  slide  chambers  containing  tissue  culture 
medium  demonstrated  that  the  energy  distribution  within  and  between  chambers  of  the 
two  chamber  model  is  nonuniform,  the  difference  between  chambers  being  about  50%.  In 
the  four  chamber  model  measurements  demonstrated  that  differences  in  energy  deposition 
remain  within  10%  of  an  average  value.  Therefore,  the  four  chamber  model  was  selected 
for  use  in  experiments.  Three  of  the  chambers  contain  tissue  culture  medium  and  cells, 
the  fourth  contains  medium  only.  A  hole  is  drilled  in  the  cover  of  the  chamber  and  a  tem¬ 
perature  probe  is  introduced.  In  this  way  each  experiment  can  be  carried  out  in  triplicate 
on  simultaneously  exposed  samples  with  simultaneous  continuous  dosimetric  control. 

Dosimetry  and  Temperature-Time  Analysis 

During  an  experiment  the  temperature  in  the  exposed  and  the  sham-exposed  samples 
is  recorded  at  regular  intervals,  with  a  minimum  interval  of  less  than  1  s.  The  temperature 
is  recorded  before  the  beginning  of,  during,  and  after  the  exposure.  The  “on”  and  “off” 
times  of  exposure  are  recorded  on  the  computer  by  the  operator  using  a  soft-key  interrupt 
capability.  At  the  conclusion  of  a  run  the  T(f),  i.e.,  temperature  (T)  versus  time  (t) 
behavior  is  analyzed  to  determine  SAR. 

If  the  SAR  is  5  (Wkg-1 )  and  the  specific  heat  of  the  sample  is  C(Jkg-1K_1 ),  then  the 
rate  of  heating  during  microwave  exposure,  for  small  temperature  variations  over  which  C 
can  be  assumed  to  be  constant,  is 


If,  in  the  absence  of  deliberate  heating,  a  sample  is  above  or  below  its  equilibrium  temper¬ 
ature,  Teg,  with  its  surroundings,  then  for  small  temperature  differences  from  equilibrium, 
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its  natural  rate  of  temperature  change,  from  Newton's  law  of  cooling,  is  of  the  form 

(f)=<r„-r),  (2) 

where  a  is  a  constant  that  depends  on  sample  geometry,  insulation,  etc.  In  general,  when 
microwave  power  is  applied  to  a  sample  that  was  not  in  thermal  equilibrium  with  its 
surroundings  at  the  start  of  exposure: 

f  -§  +  *<Th-T)  (3) 

The  solution  to  eqn.  (3)  can  be  written  in  the  form 

r-r„=(^  +  r„-7;)(i -*-»■)  (4) 

where  T0  is  the  control  temperature  at  t  =  0.  Either  eqn.  (3)  or  eqn.  (4)  can  be  used  to 
determine  the  SAR.  From  eqn.  (3), 


where  The  subscript  5  indicates  the  application  of  microwave  power  corresponding  to 
SAR  S.  Therefore,  determination  of  the  rate  of  temperature  change  before  and  during 
application  of  microwave  power,  or  during  and  after  the  application  of  microwave  power, 
determines  the  SAR. 

If  the  sample  is  in  thermal  equilibrium  before  exposure  begins,  eqn.  (4)  reduces  to 
the  simpler  form 

T  -  T.  =  A(i  _  e— )  (6) 

which  can  be  used  to  find  S  from  T(t).  In  general,  the  approach  using  eqn.  (5)  is  preferable, 
since  for  times  near  t  =  0,  the  exponential  can  be  expanded  in  quadratic  form. 

In  practice,  this  analysis  is  performed  automatically.  At  the  end  of  a  run,  the  point 
of  inflection,  or  turning  point,  in  T(t)  is  found  numerically.  This  is  done  by  checking  the 
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change  in  slope  of  the  temperature  as  a  function  of  time.  In  practice,  this  change  in  slope 
occurs  immediately  upon  the  application  or  removal  of  microwave  power,  within  the  time 
resolution  of  the  data  acquisition  system.  The  temperature  profile  to  the  left  and  right 
of  this  point  is  fitted  by  least-squares  to  a  linear  or  quadratic  function  and  the  change  in 
•lope  at  the  beginning  of  the  exposure  period  yields  the  SAR. 

In  practice  temperatures  in  the  exposed  sample  are  recorded  each  1,  2  or  3  seconds, 
the  readings  from  the  sham-exposed  sample  being  taken  at  10  times  longer  intervals.  In 
experiments  in  which  the  sample  is  microwave  heated  for  longer  periods  at  a  predetermined 
elevated  temperature,  measurements  are  made  less  frequently  once  equilibrium  is  reached, 
and  serve  only  as  control  of  the  stability  of  exposure  conditions. 

Measurements  made  in  air  or  water  in  the  incubator  outside  the  waveguides,  demon¬ 
strated  oscillations  of  ±0.2°C  resulting  from  the  on-off  cycle  of  heating  controlled  by  the 
thermostat.  Over  a  24  hour  period  such  oscillations  may  reach  ±0.5°C.  Measurements 
in  medium  in  chamber-slides  placed  within  the  waveguides  demonstrated  that  these  os¬ 
cillations  are  attenuated,  the  metal  structure  of  the  waveguide  acting  as  a  heat  sink  or 
source.  Over  a  24  hour  period  the  oscillations  do  not  exceed  0.05°C,  and  usually  remain 
within  0.02°C.  The  sensitivity  of  the  temperature  probes  is  nominally  0.01°C,  however 
comparison  of  successive  readings  permits  extrapolation  to  0.005°C.  Based  on  data  in  the 
literature  the  biological  endpoints  examined  are  not  affected  by  temperature  excursions 
within  0.5°C. 

The  software  written  for  reading  temperatures  from  two  Narda  probes  and  for  “smooth¬ 
ing  and  plotting  of  experimental  data”,  s.e.,  for  SAR  computation  is  enclosed  as  attache- 
ment  1.  The  software  was  slightly  modified  to  accommodate  the  alternate  possibility  for 
taking  reading  from  Vitek  probes. 

The  dosimetric  system  described  above  can  be  used  in  experiments  on  the  effects  of 
heating  caused  not  only  by  microwaves  but  also  by  any  means  (eg.  conventional,  ultra¬ 
sound).  Because  of  its  wide  applicability  the  system  is  the  subject  of  two  brief  communi¬ 
cations:  one  at  the  IRPA-7  Congress  in  April,  1988  in  Sydney  Australia  (attachment  2), 
and  a  second  one  at  the  10th  Annual  Meeting  of  the  BEMS  in  June,  1988  in  Stamford, 
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CT  (attachment  3).  Finally,  a  detailed  description  was  published  in  Health  Physics  56, 
303-307, 1989,  (attachment  4). 

STUDIES  OF  LYMPHOCYTE  CULTURES 
Collection  of  Blood  Samples  and  Separation  of  Lymphocytes 

Approval  for  the  use  of  human  blood  was  obtained  November  19, 1986  from  the  Human 
Subjects  Research  Review  Board,  Office  of  the  Surgeon  General,  Department  of  the  Army. 
Volunteers  were  informed  about  the  nature  of  the  study,  potential  risks  and  the  potential 
benefit  derived  from  the  study,  and  signed  the  informed  consent  form,  as  approved  by  the 
HSRRR.  20ml  of  blood  was  withdrawn  by  venipuncture  under  aseptic  conditions  into  two 
heparinized  (143  USP  units)  sterile  evacuated  blood  collection  tubes  (Vacutainer  No.  6480, 
Becton  and  Dickinson).  The  blood  was  transported  to  the  laboratory  at  room  temperature. 
Within  one  hour  from  withdrawal  mononuclear  cells  (lymphocytes  and  monocytes)  were 
separated  under  sterile  conditions  from  the  whole  blood  by  gradient  centrifugation  in 
Lymphocyte  Separation  Medium  (Boehringer,  Mannheim,  West  Germany)  according  to  the 
manufacturer’s  instructions.  Following  separation  the  cell  suspension  was  washed  twice  in 
RPMI 1640  cell  culture  medium  (Gibco).  An  aliquot  was  withdrawn  for  the  determination 
of  the  cell  concentration,  cell  viability  and  cytological  (cytocentrifuge)  preparations,  the 
remainder  was  used  for  establishing  cultures. 

The  cells  were  counted  in  a  bright  line  hemacytometer  (Neubauer  type,  Spencer).  Vi¬ 
ability  was  tested  by  the  Trypan  blue  exclusion  test,  neutral  red  stain  for  the  vacuome,  and 
Janus  green  stain  for  mitochondria.  Cytocentrifuge  (Cytospin  2  cytocentrifuge,  Shandon) 
preparations  were  air  dried,  fixed  in  methanol,  stained  by  a  combined  Wright-Giemsa  stain 
and  cell  morphology  was  analyzed  under  a  microscope  (see  below). 

The  usual  yield  of  cells  ranged  from  lxlO7  to  3xl07.  The  suspension  consisted  of 
90-95%  lymphocytes,  the  remainder  of  cells  being  monocytes  and  occasional  (below  1%) 
granulocytes.  The  viability  of  cells  ranged  from  96.5  to  99.5%,  and  was  usually  about  98%. 

Cell  Cultures 

Four-chamber  culture  chamber  slides  were  filled  under  sterile  conditions  with  1  ml  of 
chromosome  medium  1A  (Gibco)  with  or  without  phytohemagglutinin  (PHA)  depending 
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on  the  protocol  appropriate  for  the  particular  experiment.  10®  cells  were  added  to  3  of  the 
chambers,  the  4th  contained  medium  only.  A  0.2  mm  hole  was  drilled  over  the  chamber 
that  does  not  contain  cells  for  the  introduction  of  the  temperature  probe.  (For  experiments 
on  spontaneous  transformation  a  set  of  experimental  samples  consists  of  four  slides:  one 
with  PHA  (5  or  10  mg/ml)  positive  control  incubated  in  a  conventional  incubator,  and 
three  without  PHA.  Out  of  these  one  is  incubated  in  a  conventional  incubator  at  37°C, 
one  is  microwave-exposed  and  one  is  sham  exposed). 

Cell  Harvest 

Following  incubation  the  cultures  without  PHA  were  mixed  within  the  culture  cham¬ 
ber  using  a  Pasteur  pipette,  and  transferred  to  microfuge  tubes.  The  cell  count  was 
established,  viability  was  tested  and  cytocentrifuge  preparation  were  made  (as  described 

In  the  case  of  cultures  incubated  with  PHA,  lOpl  of  colcemide  solution  (10  mg/ml, 
Gibco)  were  added  to  each  chamber  2  hours  before  harvest.  At  harvest  the  cultures  were 
transferred  to  microfuge  tubes,  spun  down  at  103  rpm  for  8  min  in  an  Eppendorf  micro  cen¬ 
trifuge,  and  the  supernatant  was  discarded.  1  ml  1%  sodium  cirtrate  solution  prewarmed 
to  37°  C  was  added  with  vortexing.  The  cells  were  kept  in  this  solution  for  8  min,  spun 
down  under  the  same  conditions  as  above,  and  the  supernatant  was  discarded.  Fixative 
(cold  glacial  acetic  acid:  methanol,  1:3  vol/vol)  was  added  dropwise  with  vortexing.  After 
at  least  20  min  the  fixative  was  changed  once  or  twice.  A  droplet  of  the  cell  solution  in 
fixative  was  placed  on  a  microscope  slide  wet  with  20%  ethanol,  and  the  slides  were  flamed. 
Following  this  the  slides  were  air  dried,  and  stained  with  Giemsa  solution  in  phosphate 
buffer  pH  6.8  for  8  min  or  treated  with  trypsin  and  stained  with  Giemsa  according  to  the 
Seabright  procedure  (G-banding). 

Analysis  of  Microscope  Preparations 

Cytocentrifuge  slides  were  examined  under  a  100  x  oil  immersion  objective  and  the 
cell  and  nucleus  areas  were  determined  using  an  “Optomax”  image  analyzer.  The  original 
software  supplied  by  the  manufacturer  was  rewritten.  The  slides  was  coded  and  the  person 
who  made  the  measurements  was  not  aware  of  the  treatment  (exposed,  sham-exposed  or 
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conventional  incubator).  The  code  number  of  the  slide,  the  area  of  the  cell,  the  area  of  the 
nucleus,  the  nucleus-cytoplasmic  ratio  and  a  symbol  for  morphological  features  of  the  cell 
(S  for  a  small  lymphocyte,  1  for  an  intermediate  form  and  B  for  a  lymphoblastoid  cell) 
were  recorded,  and  stored  on  a  floppy  disk. 

The  process  of  lymphoblastoid  transformation  was  initially  evaluated  by  most  re¬ 
searchers  on  the  basis  of  cell  morphology.  Following  stimulation  the  lymphocyte  under¬ 
goes  gradual  enlargement,  accompanied  by  structural  alterations  of  the  nuclear  chromatin, 
and  becomes  an  intermediate  cell.  The  transition  of  the  intermediate  into  a  lymphoblas¬ 
toid  cell  consists  in  further  enlargement  of  the  cell,  accompanied  by  additional  changes  in 
chromatin  structure  and  the  appearance  of  nuclei.  Following  a  phase  of  DNA  synthesis 
(S-phase)  and  pre-mitotic  period  during  which  no  morphological  changes  occur  (g2  phase), 
the  lymphoblastoids  undergo  cell  division.  Evaluation  of  the  transformation  process  based 
on  qualitative  assessment  of  cell  morphology  is  subjective,  and  varies  with  different  ob¬ 
servers.  Because  of  this,  determination  of  tritiated  thymidine  incorporation  into  DNA  is  the 
most  widely  used  method  for  a  quantitative  assessment  of  lymphoblastoid  transformation. 
However,  this  method  does  not  fully  reflect  the  process  of  lymphoblastoid  transformation, 
being  only  a  measure  of  the  rate  of  DNA  synthesis.  Small  and  intermediate  lymphocytes, 
as  well  as  lymphoblastoids  in  other  stages  of  the  cell  cycle  rather  than  the  S-phase  are  not 
taken  into  account.  However,  using  this  method  it  is  not  possible  to  distinguish  between 
an  increase  or  a  decrease  in  tritiated  thymidine  incorporation  caused  by  a  change  in  the 
number  of  cells  synthesizing  DNA  or  by  a  change  in  the  rate  of  synthesis  alone.  As  the 
present  study  requires  incubation  of  cells  at  elevated  temperatures,  which  may  affect  not 
only  the  basal  metabolic  rate  but  also  the  DNA  synthesis  rate,  a  morphological  approach 
supplemented  by  objective  measurements  of  the  cell  and  nucleus  area  was  adopted. 

To  validate  this  approach,  three  different  observers  analyzed  cytocentrifuge  prepara¬ 
tions  from  cultures  incubated  with  PHA  at  three  concentrations  of  5,  7.5,  and  10  mg/m J. 
The  last  one  was  indicated  by  the  manufacturer  as  the  optimal  concentration.  The  distri¬ 
bution  of  the  size  (area)  of  the  cells  demonstrated  a  dependence  on  PHA  concentration. 
The  area  of  the  nucleus  was  distributed  randomly  and  did  not  correlate  with  the  morpho¬ 
logical  classification  of  the  cell.  Consequently  the  nucleus-cytoplasmic  ratio  could  not  be 
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used  for  the  evaluation  of  the  transformation  process. 

Cells  of  size  between  80  and  120  pi x?  were  classified  by  different  observers  as  small  or 
intermediate  lymphocytes,  and  cells  between  190  and  210  pm2  as  intermediate  or  blastoid 
forms.  Therefore  based  on  morphology  and  measurements  of  the  cell  area,  cells  up  to 
100/im3  were  grouped  in  one  class  of  small  lymphocytes.  Cells  larger  than  200/rni7  were 
considered  as  blastoid  forms.  Intermediate  lymphocytes  were  additionally  subdivided  into 
three  groups:  larger  than  100/im2  and  smaller  or  equal  to  120/im2,  larger  than  120/im2 
and  smaller  or  equal  to  150^m2,  and  larger  than  150/im2  and  smaller  or  equal  to  200 /im2. 

The  distribution  of  the  cell  size  in  preparations  from  cultures  incubated  72h  with 
PHA  at  the  three  tested  concentration  is  shown  in  Table  1.  Significant  differences  were 
found  in  the  number  of  small  lymphocytes  (area  <  100pm2)  and  lymphobl&stoids  (area  > 
200 pm2).  The  decrease  in  the  number  of  small  lymphocytes  and  the  increase  in  the  number 
of  blastoid  cells  with  increasing  PHA  concentration  validates  the  use  of  measurements  of 
the  cell  area  for  the  assessment  of  lymphocyte  transformation. 

Table  1 

The  Distribution  of  cell  size  (area,  pm2)  following  72h  incubation 
with  PHA  at  Various  Concentrations  (mg/m 1) 


PHA  mg/ml 
Area  (a)  pm2 

5.0 

wm 

10.0 

a  <100 

42.5  ±  3.7 

26.8  ±  4.5 

100  <  a  <  120 

13.4  ±  2.0 

13.0  db  2.8 

120  <  a  <  150 

15.6  ±  3.4 

12.1  ±  3.1 

10.0  ±  2.1 

150  <  a  <  200 

19.0  ±  2.6 

24.3  ±  2.8 

22.6  ±  3.0 

a  >  200 

10.6  ±  4.8 

25  ±  3.5 

46.2  ±  4.5 

Chromosomal  preparations  were  examined  using  a  100  x  oil  immersion  objective. 
Metaphase  plates  were  analyzed  for  numerical  structural  aberrations. 

Results  were  analyzed  using  tb*  “Statgraphics”  program  by  Statistical  Graphics  Cor¬ 
poration.  Analysis  of  variance  (ANOVA)  was  used  to  compare  the  numbers  of  cells  in 
particular  classes,  distributions  were  compared  using  multiple  regression  methods.  Results 
were  considered  significant  at  the  p  <  0.05  level. 
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Effects  Of  Exposure  On  Spontaneous  Transformation 

Blood  was  obtained  from  healthy  donors  by  venipuncture,  lymphocytes  were  separated 
as  described  above,  and  transferred  to  Chromosome  Medium  1A  without  PHA  (Gibco). 
The  cell  concentration  was  adjusted  to  106/ml.  One  milliliter  of  cell  suspension  was  in¬ 
troduced  into  each  of  the  3  chambers  of  the  chamber-slide,  the  4th  chamber  contained 
medium  only,  and  contained  the  sensor  for  the  dosimetric  system.  Thus  each  experiment 
was  carried  out  in  triplicate,  and  usually  repeated  3  times  (a  total  of  9  samples),  some 
experiments  were  carried  out  twice  (6  samples). 

Each  experiment  comprised  three  chamber-slides,  one  exposed,  the  second  one  sham 
exposed,  and  the  third  one  incubated  in  a  separate  incubator  set  at  37°C,  serving  as 
control. 

The  experiments  were  carried  out  in  three  series: 

•  exposed  to  CW  2450  MHz  at  a  non-heating  level,  i.e.  at  an  average  SAR  of  about  1 
W/kg  (  ranging  from  0.8  to  1.3  W/kg),  and  at  heating  levels  resulting  in  an  increase  in 
temperature  of  the  sample  by  0.5°C  (SAR  1.8-2.3  W/kg),  1  °C  (SAR  3.5-4.5  W/kg),  1.5°C 
(SAR  6. 8-8.3  W/kg),  and  2°C  (9.8-12.3  W/kg);  the  temperature  of  the  sample  varied 
over  the  incubation  period  by  ±0.2°C;  the  range  of  SARs  is  indicated  taking  into  account 
the  experimentally  observed  variations  and  the  ±10%  non- uniformity  of  energy  absorption 
over  the  three  chambers  of  the  chamber  slide; 

•  exposed  to  PW  2450  MHz  at  the  same  non-heating  and  heating  levels  of  temperature  and 
SARs  (for  these  PW  exposures  the  microwave  pulse  width  was  Ifis  and  the  pulse  repetition 
frequency  varied  from  100  to  1000  pps  as  the  exposure  changed  from  non- heating  to  the 
highest  heating  level); 

•  incubated  at  the  temperatures  obtained  by  microwave  heating  (37.5°C,  38° C,  38.5°C, 
and  39° C)  in  an  incubator  with  adjusted  thermostat,  samples  incubated  at  30° C  serving 
as  controls  and  comparison  to  the  non-heating  microwave  level. 

There  were  no  significant  differences  between  sham-exposed  and  control  cultures,  and 
therefore  the  results  from  these  were  pooled.  Table  2  presents  the  distribution  of  cell  sizes 
in  these  cultures. 
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Table  2 

The  distribution  of  cell  sizes  in  control  and 
sham-exposed  cultures  incubated  for  120h 


Area(a) 

pm2 

Number  of  cells 

meanis.d. 

a  <100 

53.7  ±  2.9 

100<  a  <120 

30.9  ±  1.8 

120<  a  <150 

10.1  ±  1.0 

150<  a  <200 

2.1  ±  0.8 

a  >200 

1.2  ±  0.9 

Table  3  presents  the  distribution  of  cell  sizes  in  conventionally  heated  cultures.  There 
are  no  significant  differences  between  cultures  incubated  at  30° C  and  37.5oC.  At  38°C  the 
number  of  small  lymphocytes  (axea<  100/jm2decreases,  and  the  number  of  cells  with  areas 
between  120^m2and  200/zm2increases,  thus  indicating  an  enhancement  of  the  transfoma- 
tion  process.  At  38.5°C  the  number  of  cells  with  an  area  between  100  and  120/im2decreases, 
and  the  number  of  cells  with  area  between  150  and  200/im2increases.  No  blastoid  cells  were 
found.  This  may  be  interpreted  as  enhancement  of  transformation  (blastoid  cells  divided 
and  the  daughter  cells  are  smaller)  or  an  inhibition.  The  latter  interpretation  is  less  likely 
because  of  the  results  obtained  at  38°C  ,  and  the  lack  of  cell  destruction  in  the  preparations. 
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Table  3.  The  distribution  of  cell  sizes  in  control  cultures 
incubated  at  37.5,  38,  38.5,  and  39°C  , 
with  conventional  heating 


Number  of  cells 

mean±  s.d. 

Incubation 

temperature 
Area  (a)  pm2 

37°  C 

37.5°C 

38°C 

38.5°C 

39°  C 

a<100 

55.2±3.8 

54.8±2.2 

39.6±3.0 

53.0±4.5 

88.7±2.2 

100<  a  <120 

32.4±2.9 

29.4±3.6 

29.6±7.1 

4.7±2.3 

9.1±1.9 

120<  a  <150 

8.8±2.0 

11.2±1.7 

23.2±3.4 

16.2±2.2 

1.6±1.4 

150<  a  <200 

1.6±0.9 

1.3±1.0 

3.4±1.0 

15.8±2.0 

0.3±0.5 

a  >200 

1.0±0.9 

2.4±1.7 

1.4±1.0 

0 

0.1±0.3 

In  all  experiments  described  above  the  total  count  of  cells  did  not  change  significantly, 
being  of  the  order  of  106/ml  with  variations  witliin  the  error  of  the  method.  Cell  viability 
also  did  not  change,  remaining  at  the  initial  level  of  94  to  98%.  At  39°C  the  number  of 
small  lymphocytes  increased  significantly,  while  the  number  of  lymphoblastoid  cells  was 
smaller  than  in  cultures  incubated  at  37,  37.5,  and  38° C  ,  the  difference,  however,  was 
not  statistically  significant.  Cellular  debris  was  present  in  the  preparations.  The  total 
cell  count  decreased  from  6.7  x  105  to  9.1  x  105  cells/ml.  Viability  of  cells  decreased,  and 
ranged  from  84  to  93%.  These  results  may  be  interpreted  as  indicative  of  heat  cell  killing. 

Table  4  presents  the  results  obtained  following  CW  2450  MHz  microwave  heating. 
No  statistically  significant  differences  were  noted  between  this  and  the  preceding  series, 
with  the  possible  exception  of  cultures  heated  to  39° C  .  At  this  temperature  signs  of  cell 
destruction  were  noted.  The  total  cell  count  ranged  from  6.4  x  105  to  8.5  x  105  cells/ml. 
Viability  varied  from  75.5  to  81.5%.  In  view  of  this,  no  reliable  determination  of  the 
distribution  of  cell  sizes  could  be  made. 
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Table  4.  The  distribution  of  cell  sizes  in  cultures  exposed  to 
CW  2450  MHz  microwaves  at  non-heating  (3T°C  )  and  heating  levels 


Number  of  cells 

mean  ±  s.d 

Incubation 

temperature 
Area  (a)  fi m2 

37°C 

37.5°C 

38°C 

38.5°C 

39°  C 

a<100 

53.3±2.1 

40.0±3.4 

39.4±3.0 

43.6±4.3 

partial  cell 

100<  a  <120 

31.7±3.0 

36.0±1.9 

30.0±4.9 

15.2±2.7 

destruction, 

120<  a  <150 

10.2±2.4 

12.2±2.0 

22.8±2.8 

17.6±1.7 

measurements 

150<  a  <200 

1.7±0.7 

4.8±1.6 

2.1±0.9 

22.6±2.8 

unreliable 

a  >200 

1.0±0.9 

2.4±1.7 

5.7±1.7 

1.0±0.9 

Table  5  presents  the  results  obtained  following  PW  2450  MHz  microwave  heating. 
Exposure  at  the  non-heating  level  results  in  an  increase  in  the  number  of  intermediate 
forms  with  an  area  greater  than  150/jm2and  smaller  or  equal  to  200/im2,  the  number  of 
blastoid  forms  increases  too.  These  increases  may  be  considered  real  in  view  of  the  results 
obtained  at  37.5  and  38°C  .  At  these  temperatures  a  statistically  significant  increase  in  the 
numbers  of  lymphoblastoid  cells  (area<200^m2)  occurs.  At  38.5°C  cell  destruction  renders 
the  determination  of  cell  size  unreliable,  and  no  measurements  were  made.  In  view  of  this, 
planned  experiments  at  39°C  were  not  performed,  no  reliable  results  were  expected. 
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Table  5.  The  distribution  of  cell  sizes  in  cultures  exposed  to 
PW  2450  MHz  microwaves  at  non-heating  (37° C  )  and  heating  levels 


Number  of  cells 

mean  ±  s.d 

Incubation 

temperature 
Area  (a)  fxm2 

37°  C 

37.5°C 

386C 

38.5°C 

39°  C 

a<100 

59.7±3.9 

44.6±3.8 

46.8±4.7 

partial 

100<  o  <120 

15.4±1.8 

20.6±2.4 

16.3±2.0 

cell 

120<  a  <150 

13.7±1.4 

16.6±1.5 

16.3±2.0 

destruction 

not 

150<  a  <200 

4.9±1.8 

3.9±1.0 

10.6±2.0 

measurements 

performed 

a  >200 

6.3±1.7 

14.4±1.2 

15.9±3.7 

unreliable 

Effects  of  Exposure  on  PHA-induced  Transformation 

Blood  was  obtained  from  healthy  donors  by  venipuncture,  lymphocytes  were  sepa¬ 
rated  as  described  above,  and  transferred  to  Chromosome  Medium  1A  containing  5  mg/ml 
PHA  (Gibco).  The  cell  concentration  was  adjusted  to  106  cells/ml.  One  milliliter  of  cell 
suspension  was  introduced  into  each  of  3  chambers  of  the  chamber  slide,  the  4th  cham¬ 
ber  contained  medium  only,  and  housed  the  sensor  of  the  dosimetric  system.  Thus  each 
experiment  was  performed  in  triplicate,  and  usually  was  repeated  three  times:  in  some 
instances  twice.  Each  experiment  comprised  three  chamber-slides,  one  exposed,  a  second 
one  sham-exposed,  and  a  third  one  incubated  at  30°C  in  a  separate  incubator,  serving  as 
control.  The  conditions  chosen  for  these  experiments  were  based  on  previously  obtained 
results.  No  experiments  were  carried  out  under  conditions  that  did  not  demonstrate  signif¬ 
icant  differences  between  experimental  and  control  samples,  or  resulted  in  cell  destruction 
rendering  measurements  unreliable.  The  experiments  were  carried  out  in  three  series: 

•  exposed  to  CW  2450  MHz  microwaves  over  the  whole  72h  incubation  period  at  a  non¬ 
heating  level,  i.e.  at  an  average  SAR  of  about  lW/kg  (ranging  from  0.7  to  1.3  W/kg) 
and  at  a  level  resulting  in  an  increase  in  the  temperature  of  the  sample  by  1°C  (38°C  ) 
i.e.  at  an  average  SAR  of  about  4 W/kg;  the  temperature  of  the  sample  varied  over  the 
course  of  the  experiment  by  ±  0.2°C  ;  the  indicated  range  of  SAR  takes  into  account  the 
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experimentally  observed  variations  in  different  chambers  of  the  chamber-slide; 

•  exposed  over  the  whole  72h  of  the  incubation  period  to  PW  2450  MHz  microwaves  at 
the  same  levels  of  average  SAR  and  temperature  as  in  the  preceding  series;  an  additional 
set  of  experiments  consisted  in  exposure  to  PW  2450  MHz  at  an  average  SAR  resulting  in 
an  increase  in  the  temperature  of  the  sample  by  1°C  over  the  first,  second  and  third  24h 
of  the  72h  incubation  period,  following  which  the  cultures  were  harvested. 

•  cultures  incubated  for  72h  in  an  incubator  adjusted  to  38°C  and  compared  to  cultures 
incubated  at  37°C  . 

The  results  of  these  experiments  are  summarized  in  Table  6. 


17 


Table  6.  The  distribution  of  cell  sizes  following  72h 
exposure  to  PW  2450  MHz  microwaves  at  a  non-heating  (37°C  ) 
and  heating  level  (38°C  )  compared  to  sham-exposed  samples. 


Number  of  cells 

mean  ±  s.d 

Incubation 

Shams  (37°C  ) 

37°  C 

38°C 

temperature 

Area  (a)  /jm2 

a<100 

43.2±7.4 

14.3±5.1 

11.7±5.3 

100<  a  <120 

14.3±4.6 

16.7±3.0 

5.9±2.9 

120<  a  <150 

14.9±5.3 

20.0±6.2 

7.8±2.7 

150<  a  <200 

11.7±5.1 

19.9±4.6 

14.5±4.1 

a  >200 

15.9±6.4 

29.1±6.6 

57.7±10.4 

Control  and  sham-exposed  cultures  did  not  differ  significantly  from  each  other.  Heat¬ 
ing  to  38°  C  by  conventional  means  and  by  CW  microwave  exposure  enhanced  transfor¬ 
mation  to  a  limited  extent.  Exposure  to  CW  microwaves  at  a  non- heating  level  did  not 
affect  the  transformation  process.  Exposure  to  PW  2450  MHz  during  the  whole  incubation 
period  of  72h  enhanced  transformation  at  both  the  heating  and  non-heating  levels.  Expo¬ 
sure  to  PW  2450  MHz  at  the  heating  (38°  C  )  level  during  the  first  24h  of  the  incubation 
period  did  not  affect  transformation,  no  significant  differences  between  exposed  and  sham- 
exposed  samples  were  demonstrated.  Exposure  during  the  second  24h  of  the  incubation 
period  resulted  in  an  inhibition  of  transformation  and  a  complete  absence  of  cell  divisions. 
Exposure  dining  the  last  24h  enhanced  transformation. 

No  chromosomal  aberrations  were  demonstrated  under  any  exposure  conditions. 

CONCLUSIONS 

Spontaneous  transformation  of  human  lymphocytes  in  vitro  is  enhanced  by  conven¬ 
tional  heating  and  CW  2450  MHz  microwave  heating.  No  significant  differences  were 
demonstrated  between  preparations  obtained  from  cultures  heated  to  the  same  tempera¬ 
ture  by  conventional  means  and  by  CW  microwaves.  PW  microwave  exposure  enhances 
transformation  at  a  non-heating  level.  Further  enhancement  of  transformation  was  seen 
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in  preparations  heated  by  PW  exposure  to  37.5  and  38° C  ,  the  number  of  transforming 
cells  being  significantly  higher  than  in  cultures  heated  by  CW  exposure.  Following  PW 
exposure  at  the  heating  level  resulting  in  an  increase  in  the  termperature  of  the  sample 
by  1.5°C  (to  38.5°C  )  or  2°C  (to  39°C  )  numerous  cells  undergo  destruction  making  the 
interpretation  of  the  preparation  unreliable. 

Effects  of  CW  exposure  seem  to  be  related  to  heating,  t.e.  depend  upon  the  increase 
in  temperature.  Effects  obtained  following  PW  exposure  differ  from  those  seen  following 
CW  exposure,  PW  exposure  enhances  transformation  at  non-heating  levels,  while  CW 
exposure  does  not.  Thus  effects  of  CW  and  PW  2450  MHz  microwaves  are  different  and 
are  not  related  to  the  average  SAR. 

PHA-induced  transformation  is  affected  by  exposure  in  a  similar  maner,  and  leads  to 
the  same  conclusions.  Additionally  this  part  of  the  study  demonstrated  the  dependence  of 
the  effects  on  the  stage  of  the  transformation  process.  Exposure  during  initial  stages  (first 
24h  of  the  72h  incubation  period)  is  without  effect.  Exposure  during  the  second  24h  of  the 
72h  incubation  period  blocks  cell  division  and  decreases  transformation,  while  exposure 
during  the  last  24h  enhances  transformation. 

Results  of  this  study  were  presented  at  the  10th  Annual  Meeting  of  the  Bioelec¬ 
tromagnetics  Society  (2  presentations)  in  Stamford,  Connecticut  in  1988  and  at  the  7th 
International  Congress  of  the  International  Radiation  Protection  Association  (IRPA)  in 
Sydney,  Australia  (2  presentations).  Copies  of  abstracts  are  included  in  this  report.  One 
paper  describing  the  exposure  system  and  dosimetry  protocols  has  already  been  published 
in  Health  Physics  (copy  included).  Two  further  manuscripts  are  in  preparation  and  will  be 
submitted  later.  Also  included  as  an  attachment  to  this  report  is  the  computer  software 
used  for  experimental  monitoring,  temperature  history  data  acquisition,  and  dosimetry. 
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Attachment  #1 


PROGRAM  FOR  COLLECTING  DATA 

10  OPTION  BASE  1 
20  DISP  "ARE  YOU  USING  VITEK  OR  NARDA  PROBES, TYPE  V  OR  N" 

30  INPUT  PRBTYPS 

qO  PR0BC1-9.6184  $  PR0BC2-41 .459  !  Calibration  factors  for  Naraa  probe  1 

50  PR0BC3-S.9121  £  PR0BC4-42 ,38b  !  Calibration  factors  for  Narda  prone  2 

60  IF  PRBTYP$-"N"  THEN  GOTO  90  !  Probe  type  Narda  -  don't  reset  caiib.  factors 

70  PROBC1-100  6  PROBC2-0  !  Calibration  factors  for  Vitek  probe 

80  PRGBC3-100  PR0BC4-0  !  Calibration  factors  for  Vitek  probe 

90  ON  KEY#  1 ."POWER  ON"  GOSUB  1380 

100  ON  KEY#  4, "NEW  0BS  INT"  GOTO  1100 

110  ON  KEY#  7. "POWER  OFF"  GOSUB  1390 

120  MASS  STORAGE  IS  ":D500" 

130  DIM  ATEMP(IOOO) , ATIME< 1000) . BTEMPC 1000) ,BTIMEi 1 000) , PONT (50) ,POFFT<50) 

140  DIM  ANUM(20) ,BNUM(20) ,IPON(20> ,IPOFF<20) 

150  CLEAR 

160  !  BEEP  40,100  $  WAIT  50  »  BEEP  40,100  ®  WAIT  50  »  BEEP  40,100  $  WAIT  50  *  BE 
EP  80,150 

170  DISP  "ENTER  DAY .MONTH , YEAR" 

180  INPUT  EXPDATES 
190  EXPDATE-0 

200  DISP  "ENTER  SAMPLE  IDENTIFIER" 

210  INPUT  SAMPLES 

220  EXPDATAS-EXPDATES&SAMPLES 

230  DISP  "ENTER  TIME  OF  DAY  -  24  HOUR  CLOCK"  %  INPUT  TODS 

240  ON  ERROR  GOTO  300 

250  CREATE  EXPDATAS6" :D501 ",  1 ,32000 

260  ASSIGN#  1  TO  EXPDATAS&" :D50 1 " 

270  !  PRINT#  1  ;  EXPDATES 
280  !  PRINT#  1  ;  SAMPLES 
290  GOTO  390 

300  IF  ERRN  -128  THEN  DISP  "DISC  IS  FULL  -  INSERT  A  NEW  DISC  AND  START  AGAIN" 

310  IF  ERRN  -128  THEN  GOTO  340 
320  IF  ERRN  -63  THEN  GOTO  350 

330  DISP  "AM  ERROR  HAS  OCCURRED  -  PRESS  CONTINUE  TO  START  AGAIN" 

340  PAUSE  6'  GOTO  150 

350  DISP  "THIS  DATA  FILE  NAME  ALREADY  EXISTS  -  ENTER  ANOTHER  SAMPLE  NUMBER" 

360  INPUT  SAMPLES 

370  EXPDATAS-EXPDATES&SAMPLES 

380  GOTO  250 

390  KK-0  !  Counter  for  new  observation  times  for  probes  A  and  B. 

400  BIN-0  !  Holds  previous  number  of  observations  from  probe  A. 

410  NDPTS-0  !  Counter  for  number  of  data  points  recorded. 

420  CLEAR 

430  DISP  "This  program  will  store  1000  temps,  and  times  for  both  probes." 

440  DISP  "ENTER  TOTAL  OBSERVATION  TIME  IN  SECONDS." 

450  INPUT  OBSTIME 

460- DISP  "Probe  A  should  be  in  the  exposed/experimental  culture." 

470  DISP 

480  DISP  "ENTER  OBSERVATION  INTERVAL  FOR  PROBE  A  IN  SECONDS." 

490- INPUT  NA 

500  IF  NA>QBSTIME  THEN  520 
510  GOTO  540 

520  DISP  "You  selected  an  obs.  interval  >  the  total  obs.  time;  try  again." 

530  GOTO  430 

540  AJQBS-INT  (OBSTIME/NA) 

550  NA-NA-1000 
560  HOLD-AJOBS+BIN 
570  IF  HQLD>1 000  THEN  590 
580  GOTO  610 

590  DISP  "You  have  requested  too  many  readings  of  probe  A;  try  again." 

600  GOTO  430 


610 
620 
630 
640 
650 
660 
670 
680 
690 
700 
710 
720 
730 
740 
750 
760 
770 
780 
790 
80  0 
810 
820 
830 
840 
850 
860 
870 


MThe 

“of. 


observation 
and  ideally 


interval  for  probe  B 
<  or  ■  10  times  that 


should  be 
for  probe 


an 

ft.' 


DISP 
DISP 
DISP 

DISP  "ENTER  OBSERVATION  INTERVAL  FOR  PROBE  B  IN  SECONDS*" 
INPUT  NB 

BJOBS-INT  (OBSTIME/NB) 

NB-NB-1000 

IF  FP  <NB/NA>#0  THEN  700 
GOTO  720 

DISP  "The  obs.  time  for 
GOTO  640 

REMOTE  502  !  DVM  to  remote. 

CLEAR  3  KEY  LABEL 
REMOTE  507  f  Relay  to  remote. 

LOCAL  LOCKOUT  5  !  Local  lockout  for  DVM  and 
OUTPUT  502  ;"F0R2T0X"  !  DVM  to  DCV,2V  range 
OUTPUT  507  "A12"  !  Connect  all  AUA2 
ENTER  502  *,  AS  f  Trial  read  of  DVM. 

IJ1*0  !  Set  counter  for  number  of  times 
IJ2*0  !  Set  counter  for  number  of  times 
J*1  !  Set  counter  for  reading  probe  B. 


integer  multiple" 


B  is  not  a  multiple  of  that  for  A;  try  again. 


relay . 

__  _  .  and  cont.  on  talk, 
terminals  to  C1&C2  terminals. 


power 

power 


turned 

turned 


on . 
off. 


only. 

arrays. 


NT-NB/NA  ?  Set  interval  between  reading  probe  B. 

IF  KK>0  THEN  880 

SETTIME  0,0  !  Zero  clock  for  easy  interpretation 
!  Loop,  read  probe  A,  read  time,  put  values  into 
NNA* 1  *  NNB= 1 

FOR  1-1  TO  AJOBS  3  ENTER  502  ;  AS®  ATIME<I>-TIME  ♦(DATE  -EXPDATE>*86400  3  A 
£MP(I)-VAL  (ASI51 > 

^TEMP( I )-PR0BCl*ATEMP(I)+PR0BC2  !  Evaluate  temp,  from  voltage  reading. 

890  NNA-NNA+1 

900  DISP  "ATEMP-" ,ATEMP(I> ,"TIME*", ATIME(I) 

910  IF  I«J*NT  THEN  930  !  Check  to  see  if  probe  B  is  to  be  read. 

920  GOTO  1010 

930  OUTPUT  507  *'B1"  !  Switch  B1  to  Cl  prior  to  reading  probe  B. 

940  WAIT  200  !  Wait  200ms  for  DVM  to  read  probe  B. 

§59  l  Probe  B,  read  time,  put  values  into  arrays.  Switch  Al  to  Cl. 

|6OUTPUTR5072"4l»$€i  BTIME(J)*TIME  -<DATE  -EXPDATE >*86400  3  BTEMP<J)-VAL  <B$I51) 

^80  NNB^NNB*  •  Evaluate  temp,  from  voltage  reading. 

990  DISP  "BTEMP-" ,BTEMP( J) , "TIME-",BTIME( J)  3  DISP  3  DISP  3  DISP  3  KEY  LABEL 


1000 
1010 
1020 
1030 
1040 
1Q50 
1060 
1070 
1Q80 
1090 
1100 
1110 
1120 
1130 
1140 
1150 
1160 
1170 
1180 
1 1SO 
1200 


JmJ+1  !  Increment  counter  on  read  probe  B. 

WAIT  NA  !  Wait  time  between  reading  probe  A. 

NEXT  I  !  Increment  loop. 

LOCAL  5 

OUTPUT  702  "B12"  !  Simply  turn  off  lights  on  relay. 

PJSP  data  on  disc.  You  will  be  prompted  for  more  inputs  shortly." 

DISP  3  DISP  3  DISP 

WAIT  2000 

GDSUB  DISCPRINT 

GOTO  1400 

CLEAR. 

StIS  m?°  n?t  press  any  keys  until  the  computer  responds  with  a  beep  and" 
nteo  displays  'Data  storage  is  complete' 

PISP  key  in  new  observation  data." 

WAIT  5000 
GOSUB  DISCPRINT 
GOTO  440 

!  Write  all  data  to  disc  for  this  observation  interval 
NNA-1 
NNA-1 

ATEMP( I ) ,ATIME< I ) 


You  will  then  be  prompted  to’ 


DISCPRINT: 
PRINT#  1  ; 
FOR  1-1  TO 
PRINT#  1  ; 


T2I0  NEXT  I 
1220  PRINT*  1 


1230  FOR  1-1  TO  NNB-1 


1240  PRINT*  1 
1250  NEXT  I 
1260  PRINT*  1 


NNB-1 


BTEMP<I),BTIME<I> 


1260  PRINT*  1  ;  IJ1  !  The  number  of  times  power  turned  on. 

127.0  IF  IJ1-0  THEN  1290 

1280  FOR  I-!  TO  IJ1  $  PRINT*  1  ;  P0NT<I>  $  NEXT  I 

1290  PRINT*  1  :  IJ2  !  The  number  of  times  power  turned  off. 

1300  IF  IJ2-0  THEN  1320 

1310  FOR  1-1  TO  IJ2  9  PRINT*  1  ;  P0FFT<I>  9  NEXT  I 
1320  PRINT*  1  ;  TODS  !  Writes  time  of  day  to  disc 
1330  KK-KK+1 
1340  BIN-AJOBS 

1350  NDPTS-NDPTS+2*(NNA+NNB)+IJ1+IJ2+10 

1360  BEEP  9  CLEAR  9  DISP  "Data  storage  is  complete.”  9  RETURN 
1370  GOTO  1400 

1380  BEEP  9  IJ1-IJ1+1  9  DISP  “Power  on”  9  P0NT< IJ1 )-TIME  9  RETURN 

1390  BEEP  9  IJ2-IJ2+1  9  DISP  “Power  off”  9  POFFTC IJ2)-TIME  9  RETURN 

1400  DISP  "Do  you  wish  to  record  more  data?  ENTER  <Y)ES  or  <N)0." 

1410  INPUT  REPLYS 

1420  IF  REPLYS-" Y”  THEN  420 

1430  ASSIGN*  1  TO  * 

1440  DISP  "The  final  disc  management  is  now  being  completed" 

1450  ASSIGN*  1  TO  EXPDATASA" :D501 " 

1460  IA1-1  9  IATOT-O  9  IB1-1  9  IBTOT-O  9  I0N1-1  9  IDNT-0  9  I0FF1-1  9  IOFFT-O 
1470  FOR  J1-1  TO  KK 

1480  READ*  1  ;  ANUM(J1)9  IAT0T-IATGT+ANUM< J1 ) 

1490  FOR  I-IA1  TO  IATOT  9  READ*  1  ;  ATEHP(I) ,ATIME(I)9  NEXT  I 
1500  READ*  1  ;  BNUM(J1)9  IBT0T-IBT0T+BNUM< J1 > 

1510  FOR  I-IB1  TO  IBT0T  9  READ*  1  ;  BTEMP(I) ,BTIME(I)9  NEXT  I 

o2U  Rtriu*  i  ;  ii-UNij|>9  iON»-iUNi>iPUN( J1  >  9  IF  IPDN(J1)-0  THEN  1540 

1530  FOR  I-I0N1  TO  ION  I  w  KtHu*  i  :  ruN7<I)<e-  NEXT  I 

1540  READ*  1  ;  IP0FF(J1)9  IGFFT«IGFFT+IPOFF< J1 )  9  IF  IPOFF(J1>-0  THEN  1560 

1550  FOR  I-I0FF1  TO  IOFFT  9  READ*  1  ;  P0FFT1D9  NEXT  I 

1560  IA1 -IA1 +IAT0T  9  1B1 -IB1 +IBTQT  9  I0N1 -I0N1 +I0NT  9  I0FF1 -IGFF1 +I0FFT 

1570  NEXT  J1 

1580  READ*  1  ;  TODS 

1590  ASSIGN*  1  TO  * 

1600  PURGE  EXPDATASA" :D501" 

1610  PACK  " :D501 " 

1620  NDPTS-8-NDPTS 

1630  CREATE  EXPDATAS&" :D501 ",1 ,NDPTS 
1640  ASSIGN*  1  TO  EXPDATASA":D50 1 " 


1650  PRINT*  1 


1660  FOR  1-1  TO  IATOT  9  PRINT*  1  ;  ATEMP(I) ,ATIME(I)  9  NEXT  I 


1670  PRINT*  1 


1680  FOR  1-1  TO  IBTOT  9  PRINT*  1  ;  BTEMP< I) ,BTIME(I>  9  NEXT  I 


1680  PRINT*  1  ;  10NT  9  DISP  "MW  POWER  TURNED  ON  ",I0NT,"  TIMES"  9  FOR  1-1  TO  101 
T  9  PRINT*  1  ;  PONT<I>  9  DISP  PQNT<I>  9  NEXT  I 

1700  PRINT*  1  ;  IOFFT  9  DISP  "MW  POWER  TURNED  OFF  ", IOFFT,"  TIMES"  9  FOR  1-1  TO 
IOFFT  9  PRINT*  1  ;  POFFTU)  9  DISP  P0FFT(I>  9  NEXT  1 
1710  PRINT*  1  ;  TODS  9  DISP  "Time  of  day'\TODS 
1720  ASSIGN*  1  TO 

1730  DISP  "Sample  number  is  ", SAMPLES 
1740  DISP  "Date  of  exposure  is  ".EXPDATES 

inished  the  disc  management.  GOODBYE,  HAVE  A  NICE  DAY!" 

■  rou  END 


IATOT  9  DISP  "TOTAL  PROBE  A  READINGS  -  ", IATOT 


IBTOT  9  DISP  "TOTAL  PROBE  B  READINGS  -  ", IBTOT 


i sc  management.  GOODBYE,  HAVE  A  NICE  DAY!’ 


PROGRAM  FOR  DATA  ANALYSIS,  CURVE  SMOOTHING  AND  GRAPHICS  version  A 


10  !  . . . . . 

20  !  *  * 

30  *  *  PROGRAM  TO  READ  TEMPERATURE  READINGS  TAKEN  FROM  * 

40  !  *  TWO  NARDA  PROBES  AND  STORED  ON  DISC,  PLOT  DATA  * 

50  !  *  IN  VARIOUS  FORMS  AND  DO  DOSIMETRY  CURVE  FITTING  » 

60  !  *  * 

70  !  *  WRITTEN  BY  CHRISTOPHER  C  DAVIS,  DECEMBER,  1986  * 

80  !  * 

90  f  •*««*»«••*•••**«*••«•«•»«•»»•««•«•««••»***•«•««»••« 

100  OPTION  BASE  1 

110  MASS  STORAGE  IS  ":D500" 

120  DIM  ATEMPU 000 >,ATIME< 1000), BTEMP<1 000), BTIME( 1000), PONT<50> ,POFFT(50) 
130  DIM  ANUM(20) ,BNUM(20) ,IPON(20> ,IPOFF (20) , DUMMY! [50] 

140  ♦  ***  FUNCTION  TO  GIVE  INVERSE  VIDEO  DISPLAYS 
150  DEF  FNIVIDS(DUMMYS) 

160  LENGTH-LEN  (DUMMYS) 

170  FOR  II-l  TO  LENGTH 

180  DUMMYS t II , IIl-CHRS  (NUM  (DUMMYSt II , II ] >+128) 

190  FNIVIDS-DUMMYS 
200  NEXT  II 
210  FN  END 
220  CLEAR 

230  DISP  FNIVIDSC’DATA  DATE'?") 

240  DISP  FNIVIDS("DATE  FORMAT") 


250  DISP  FNIVID$< "DAY  MONTH, YEAR") 

260  INPUT  EXPDATES 

270  DISP  FNIVID$("ENTER  SAMPLE  ID") 

280  INPUT  SAMPLES 
290  EXPDATAS-EXPDATES&SAMPLES 
300  ON  ERROR  GOTO  320 
310  GOTO  370 

320  IF  ERRN  -67  THEN  DISP  "YOU  HAVE  ENTERED  THE  WRONG  FILE  NAME" 

330  IF  ERRN  *67  THEN  DISP  "AN  ERROR  HAS  OCCURRED  -  START  AGAIN" 

340  IF  ERRN  *67  THEN  GOTO  10 

350  CAT  " :D501 " 

360  DISP  "THESE  ARE  THE  DATA  FILES  -  CHOOSE  THE  ONE  YOU  WANT"  @  GOTO  230 
370  ASSIGN*  1  TO  EXPDATAS&" :D501 " 

380  NUMSMO-O  ®  NTEST-Q  $  SAR-0  @  T0DS«"0" 

390  IATOT-O  &  IBTOT-O  €>  IONT-O  G  IDFFT-0  @  IA1-1  6  IB1-1  6  IDN1-1  %  I0FF1-1 
IATOT  !  READ  NO  OF  PROBE  A  READINGS 

TO  IATOT  %  READ*  1  ;  ATEMP(I) , ATIME(I)#  NEXT  I  !  PROBE  A  DATA 
IBTOT  !  READ  NO  OF  PROBE  B  READINGS 

TO  IBTOT  &  READ*  1  ;  BTEMP(I)  ,BTIME<I>€>  NEXT  I  !  PROBE  B  DATA 
.  IONT  f  NO  OF  TIMES  MW  PWR  TURNED  ON 
450  FOR  I-ION1  TO  IONT  S  READ*  1  ;  PONT<I)S  NEXT  I  !  ON  TIMES 
460  READ*  1  ;  IOFFT  !  NO  OF  TIMES  MW  PWR  TURNED  OFF 


400  READ*  1 
410  FOR  I-IA1 
420  READ*  1  ; 
430  FOR  I-IB1 
440  READ*  l 


PQFFT<I>&  NEXT  I  !  OFF  TIMES 


520 

530 

540 


',  IATOT 


US  J;?°KoJ°.I9FFT  • READ' 1  * - -  —  *  .  u.r 

IA1-IA1+IATOT  ®  IB1 -IB1 +IBTOT  »  ION1 -IQN1 +IONT  9  IOFF1 -IOFF 1 +IOFFT 
480  ON  ERROR  GOTO  499 
491  READ*  1  ;  TODS 

499  OFF  ERROR  «  ASSIGN*  1  TO  * 

500  !  PURGE  EXPDATAS4" :D501" 

510  DISP  "TOTAL  PROBE  A  READINGS  -  _ _ 

!pgR  1-1  T0  IAT0T  ®  DISP  ATEMP(I),ATIME<I),I  »  NEXT  I 
DISP  "TOTAL  PROBE  B  READINGS  -  ", IBTOT 
ten  T0  IBT0T  •  disp  BTEMP<I),BTIME<I),I  »  next  i 

trS  2JSP  Mw  POWER  TURNED  ON  ",IONT,"  TIMES" 

15?  SiSP  ‘’Mw  P0WER  BURNED  OFF  ",  IOFFT,"  TIMES" 
loS  ^IPt’7HLmw  PWR  WAS  TURNED  ON  AT  THE  FOLLOWING  TIMES" 

||J  F?R  I:.t  TO  IONT  C  DISP  PONT < I  >  $  NEXT  I 
590  DISP  Sample  ID  is  " ,FNIVIDS(SAMPLE$> 


600  DISP  "Date  of  exposure  was  " ,FNIVID$(EXPDATE$> 

601  IF  TOD$-"0”  THEN  GOTO  610 

602  DISP  "Tune  of  day  was  ",FNIVID$<TOD$> 

61 0  !  . . . . * . . . 

620  !  *  SECTION  FOR  SMOOTHING 

630  !  •  AND  PLOTTING  OF  EXPERIMENTAL  DATA 

640  !  *  * 

650  !  *  WRITTEN  BY  Christopher  C.  Davis,  June  1984.  * 

660  !  *  MODIFIED  DECEMBER  1986  * 

670’  !  *  *  ' 

680  ! . . . * . . . . 

690  ! 

700  ! 

710  DIM  AA(5) ,SN(21 ,21 > ,R <20 ) , HEAT (5) ,C0DL(5>  !  ARRAYS  USED  IN  CURVE  SMOOTHIN 
G 

720  NYSLS-"N” 

730  IF  ATIME< IATOT)>BTIME< I3T0T)  THEN  XMAX-ATIME(IATOT)  ELSE  XMAX-BTIMEdBTuT 
)  !  SET  X  AXIS  MAX  TO  MAX  OBSERVATION  TIME 
740  XMIN-0  ♦  INITIALLY  SET  X  ORIGIN  TO  0 
750  JMIN-1  6  KMIN= 1 

760  DISP  "TEMPERATURE  MEASUREMENT  INTERVAL  WAS",XMAX 

770  YMIN-10'20  «  YMAX— <10~2G> 

780  !  ***  FIND  MAXIMA  AND  MINIMA  OF  DATA 

790  FOR  1-1  TO  IATOT 

800  IF  ATEMP(I>>YMAX  THEN  YMAX-ATEMP(I) 

810  IF  ATEMPdXYMIN  THEN  YMIN-ATEMPd) 

820  NEXT  I 

830  FOR  1-1  TO  IBTOT 

840  IF  BTEMP( I)>YMAX  THEN  YMAX-BTEMP(I) 

850  IF  BTEMPdXYMIN  THEN  YMIN-BTEMPd) 

860  NEXT  I 

870  NUMPLT-0  9  IF  NTEST#Q  THEN  CLEAR 

880  DISP  "DO  YOU  WISH  TO  PLOT  PROBE  A,  PROBE  B,  OR  BOTH?  "&FN IV IDS < "ENTER  T .2 
OR  3") 

890  DISP  "YOU  CAN  ONLY  CALCULATE  SAR  USING  PROBE  A" 

900  INPUT  NPROBE 

910  IF  NPROBE- 1  THEN  GOTO  930 

920  GOTO  940 

930  FOR  1-1  TO  IATOT  9  BTIME(I)-ATIMEd)  $  BTEMP<I)-ATEMP<I>  9  NEXT  I 
940  IF  NTEST-0  THEN  GOTO  960 

950  NYSLS«"Y" 

960  DISP  "ENTER  TIME  REGION  YOU  WISH  TO  ANALYZE" 

970  INPUT  TMIN.TMAX 

980  FOR  I-t  TO  IATOT  !  FIND  RANGE  OF  TEMP  ARRAY  TO  BE  PLOTTED 
990  IF  TMIN>ATIME<I)  THEN  JMIN-I 

1000  NEXT  I 

1010  FOR  1-1  TO  IBTOT  9  IF  TMIN>BTIMEU>  THEN  KMIN-I 

1020  NEXT  I 


1030  IF  ATIME< JMINXBTIME(KMIN)  THEN  TMIN-ATIME< JMIN)  ELSE  TMIN-BTIME(KMIN) 

1040  FOR  1-1  TO  IATOT  6  IF  TMAX>ATIME<I >  THEN  JMAX-I 
1050  NEXT  I 

1060  FOR  1-1  TO  IBTOT  9  IF  TMAX>BTIME<I>  THEN  KMAX-I 
1070  NEXT  I 


1080 
1090 
1100 
1110 
1  120 
1130 
1140 
1150 


IF  ATIME< JMAX>>BTIHE(KMAX>  THEN  TMAX-ATIHE< JMAX>  ELSE  TMAX-BT IME < KMAX ) 
IF  NTEST-0  THEN  GOTO  1400 
IF  NYSL$-,,N"  OR  NYSL$«"NO“  THEN  GOTO  1400 
***  SECTION  FOR  SELECTION  OF  CURVE  REGION  FOR  SMOOTHING 
TEMAX--1 . E20  9  TEMIN-1 .E20 

FOR  I- JMIN  TO  JMAX  9  IF  ATEMP<I XTEMIN  THEN  TEMIN-ATEMPd ) 

IF  ATEMPd )-TEMIN  THEN  JTURN-I  f  Find  inflection  point  on  temp  curve 
NEXT  I 


1160  DISP  “NUMBER  OF  POINTS  ON  COOLING  SLOPE-  JTURN-JMIN+1 

1170  DISP  “NUMBER  OF  POINTS  ON  HEATING  SLOPE-  ", JMAX-JTURN+1 

1180  DISP  “IS  THIS  SATISFACTORY?  Y(ES>  OR  N(0>“  9  INPUT  NYSATS 

1190  IF  NYSAT$«“N"  OR  NYSATS-“NO”  THEN  GOTO  960 

1200  IF  NYSLS-“N"  OR  NYSL$-“NO"  THEN  GOTO  1400 

1210  CLEAR  9  DISP  “COOLING  REGION  POINTS", JMIN, JTURN 

1220  DISP  “HEATING  REGION  POINTS” ,JTURN ,JMAX 

1230  ON  KEY#  1 /‘COOLING"  GOTO  1280 

1240  ON  KEY#  2 /‘HEATING”  GOTO  1340 

1250  DISP  FNIVIDS< "SELECT  OPTION”)  9  DISP  &  DISP 

1260  KEY  LABEL 

1270  GOTO  1270 

1280  NSM-1  $  DISP  "ENTER  POINTS  RANGE  FOR  COOLING  SLOPE" 

1290  DISP  “ENTER  0,0  FOR  HELP  IN  FINDING  TURNING  POINT  IN  HEATING  CURVE"  9  INPUT 
JMIN, JTURN 

1300  IF  JMIN-0  AND  JTURN-0  THEN  GOSUB  Find 
1310  DISP  "JTURN*" , JTURN 

1320  DISP  “ENTER  JMIN, JTURN"  9  INPUT  JMIN, JTURN 
1330  GOTO  1390 

1340  NSM=2  *  DISP  “ENTER  POINTS  RANGE  FOR  HEATING  SLOPE" 

1350  DISP  “ENTER  0,0  FOR  HELP  IN  FINDING  TURNING  POINT  IN  HEATING  CURVE"  %  INPUT 
JTURN ,JMAX 

1360  IF  JTURN-0  AND  JMAX-0  THEN  GOSUB  Find 
1370  DISP  “JTURN=", JTURN 

1380  DISP  "ENTER  JTURN, JMAX"  9  INPUT  JTURN ,JMAX 

1390  NUMPLT- 1  9  GOSUB  Smooth 

1400  GOSUB  Aplot 

1410  ON  KEY#  1 ."AGAIN"  GOTO  1530 

1420  ON  KEY#  2,"D0SE“  GOTO  1580 

1430  ON  KEY#  3,"PL0T“  GOTO  1590 

1440  ON  KEY#  4,"SCRDMP“  GOTO  3910 

1450  ON  KEY#  5  /‘EXIT"  GOTO  3920 

1460  CLEAR  9  KEY  LABEL 

1470  DISP  “SELECT  OPTION  -  "AFNIVID$<"DOSE">&"  ALLOWS  CALCULATION  OF  SAR" 

1480  DISP  “  -  "&FNIVIDS("AGAIN">&“  RE-READS  DATA  FROM  DISC” 

1490  DISP  “  -  "&FNIVID$("PLOT")&"  RE-PLOTS  LAST  PLOT  -  USE  WITH  C 

ARE" 

1500  DISP  “  -  “&FNIVID5<"SCRDMP”>4“  DUMPS  GRAPHICS  TO  PRINTER" 

1510  DISP  “  -  “4FNIVIDS("EXIT">4"  TERMINATES  EXECUTION" 

1520  GOTO  1520 

1530  CLEAR 


1540  DISP  "IF  YOU  CONTINUE  DATA  WILL  BE  RE-READ  FROM  DISC" 

1550  DISP  "PRESS  "4FNIVID$(“AGAIN“)&“KEY  IF  YOU  WISH  TO  DO  THIS" 
1560  ON  KEY#  1 /‘AGAIN"  GOTO  370  9  KEY  LABEL 
1570  GOTO  1570 

1580  GOTO  930 

1590  IF  NUMSMO-O  THEN  GOTO  1620 

1600  GOSUB  Aplot  9  GOTO  1520 
1610  GOTO  1410 

1620  DISP  “OPTION  NOT  AVAILABLE"  9  GOTO  370 
1630  Smooth: 

1640  {  . . 

1650  !  «*  «* 

1660  !  *  CURVE  SMOOTHING  SECTION 

1670  !  *  « 

1680  !  . . 


1G90 

1700 

1710 

1720 

1730 


DISP  "CURVE  FITTING  IN  PROGRESS  “4FNIVIDS (“PLEASE  WAIT") 

NUMSMO-J 

IF  NSM-2  THEN  GOTO  1750 
N- JTURN-JMIN+ 1 

FOR  1-1  TO  N  9  BTIME< I)-BTIME< JMIN+I- 1 )  9  BTEMP( I)-BTEMP< JMIN+I- 1 >  9  NEX 


T  I 
1 740 
1750 
1760 
RN+I-1) 
1770 
1780 
1790 
1800 
1810 
1820 
1830 
1840 
1850 
1860 
1870 
1880 
1890 
1900 
1910 
1920 
1930 
1940 
1950 
1960 
1970 
1980 
1990 
N) 


GOTO  1780 
N-JMAX-JTURN+1 

FOR  1-1  TO  JMAX-JTURN+1  3  BTIME<I>-BTIME< JTURN+I-1 )  *  BTEMP<I >-BTEMP< JTU 


NEXT  I 

FOR  1-1  TO  5 

R<I)-0  »  AA<I)-0  3  FOR  J-1  TO 

DISP  "WHAT  ORDER  OF  POLYNOMIAL 

INPUT  NS 

FOR  J-1  TO  NS+1 

FOR  1-1  TO  N 


21  3 
DO 


SN(I,J)-0  3  NEXT  J  $  NEXT  I 
YOU  WISH  TO  FIT  -  1 -LINEAR... 


IF  J-NS+1  AND  BTIM£<I)-0  THEN  GOTO  1860 
R(J)-BT£MP(I)*BTIMEd)  A(NS+1-J)+R(J)  $  GOTO  1870 
R< J)-BTEMP<I)+R< J) 

NEXT  I 
NEXT  J 


FOR  K-1  TO  NS+1  3  FOR  J-1  TO  NS+1  3  FOR  1-1  TO  N 
IF  2*NS- J-K+2-0  AND  BTIME(I)-0  THEN  GDTD  1920 
SN(J,K)-BTIME<I> ‘<2*N$-J-K+2)+SN<J,K>  3  GOTO  1930 
SN  <  J , K  > - 1 +SN  <  J , K ) 

NEXT  I 
NEXT  J 
NEXT  K 


GDSUB  Invt 

FOR  1-1  TO  NS+1  $  FOR  J-1  TO  NS+1 

AA<NS+2-I)-SN(NS+2-I,J)*R(J)+AA<NS+2-I>  3  NEXT  J  3  NEXT  I 

IF  NSM-1  THEN  DX=ATIME< JTURN)-ATIME< JMIN)  ELSE  DX-AT IME ( JMAX ) -ATIME ( JTUR 


2000  IF  NSM-1  THEN  XSTART-ATIMEC JMIN)  ELSE  XSTART-ATIME( JTURN) 

2010  FOR  K-1  TO  200  *  I-K+800  3  BTIME<I>-<K-1 >*DX/200+XSTART  3  BTEMP(I)-0  3  F 

OR  J-1  TO  NS+1 

2020  IF  NS+1-J-0  AND  BTIME<I)-0  THEN  GOTO  2050 

2030  BTEMP(I)-AA< J)*BTIME<I>~<NS+1-J>+BTEMP<I> 

2040  GOTO  2060 

2050  BTEMP<I)-AA( J)+BTEMP<I) 

2060  NEXT  J 

2070  NEXT  K 

2080  RETURN 

2090  Invt: 

2100  ! 

2110  ♦  MATRIX  INVERSION  BY  GAUSS-JORDAN  ELIMINATION  WITH  FULL  PIVOTING 

2120  !  SQUARE  MATRIX  OF  ORDER  NP 

2130  ! 

2140  DIM  KP(5> ,  JP(5) 

2150  NP-NS+1 

21-60  NCOL-NP 

2170  ! 

2180  !  INITIALIZE  PERMUTATION  VECTORS 

21-90  ! 

2200  FOR  1-1  TO  NP  3  KP<I)-I  3  JP<I)-I  3  NEXT  I 
2210  1 

2220  !  INVERT  MATRIX 

2230  ! 

2240  FOR  IR-J  TO  NP 
2250  ! 

2260  !  FIND  THE  PIVOT  ELEMENT 

2270  ! 


2280 

2290 

2300 

2310 


PIVOT-O 

FOR  I-IR  TO  NP  3  FDR  J-IR  TO  NP  3  K-KP(I)  3  L-JP(J) 
TEST-ABS  <SN<K,L))-ABS  (PIVOT) 

IF  TEST<-  0  THEN  GOTO  2370 


2320  IF  TEST>0  THEN  GOTO  2340 

2330  GOTO  2370 

2340  IK- I 

2350  JK-J 

2360  PIVOT-SN<K ,L) 

2370  NEXT  J  $  NEXT  I 
2380  ! 

2390  !  THE  PIVOT  ELEMENT  IS  SN(K,L> 

2400  ! 

2410  !  UPDATE  PIVOT  VECTORS 

2420  f 

2430  K-KP< IK)  £  KP(IK)-KP(IR)  0  KP(1R>-K  0  L-JP(JK)  #  JP< JK)-JP(IR)  0  JP(IR) 
L 

2440  ! 

2450  !  COMPUTE  NEW  ELEMENTS  OF  PIVOTAL  ROW 

2460  ! 

2470  FOR  J-t  TO  NCOL  0  SN(K, J)-SN(K, J)/PIVOT  0  NEXT  J 

2480  ! 

2490  !  COMPUTE  REMAINING  ELEMENTS  OF  IR'TH  STEP 

2500  ! 

2510  FOR  1-1  TO  NP 

2520  IF  I-K-0  THEN  GOTO  2590 

2530  AIL-SN( I ,L) 

2540  SN(I.L)— (AIL/PIVOT) 

2550  FOR  J-1  TO  NCOL 

2560  IF  J-L-0  THEN  GOTO  2580 

2570  SN<I,J)-SN(I,J)-AIL*SM(K,J) 

2580  NEXT  J 

2590  NEXT  I 

2600  ! 

2610  ! 

2620  !  INVERT  PIVOT  ELEMENT 

2630  ! 

2640  SN<K,L)-1 /PIVOT  0  NEXT  IR 

2650  RETURN 

2660  Aplot: 

2670  !  ******•*••******.**...«.** 

2680  !  « 

2690  !  *  PLOTTING  SUBROUTINE  « 

2700  !  **.***»*****«************ 

2710  f 

2720  NLT-1 


2730  CLEAR  0  IF  NUMPLT-0  THEN  GOTO  2930 
2740  ALPHA  0  DISP  "PARAMETERS  OF  FIT"  0  DISP 
2750  FOR  1-1  TO  NP  0  DISP  AA(I)  0  NEXT  I 
2/60  IF  NSM-2  THEN  GOTO  2790 
2770  FOR  1-1  TO  NP  0  COOL( I)-AA< I)  0  NEXT  I 
2780  GOTO  2800 

2790  FOR  1-1  TO  NP  0  HEAT(I)-AA(I)  0  NEXT  I 
2800  CHI2-0 

2810  IF  NTE$T<2  THEN  GOTO  2850 


2820  SAR-4 ,2*(HEAT ( 1 )-COOL< 1 >>-1000 

2830  IF  NS-2  THEN  SAR-2*SAR-ATIME( JTURN>+4 ,2*(HEAT<2)-C00L<2> > 
2840  DISP  FNIVIDSCSAR  IS")  , SAR , "W/kg" 

2850  FOR  1-1  TO  N  0  Y2-0  0  FOR  J-1  TO  NS*1 
2860  IF  NS+1-J-0  AND  BTIME<I)-0  THEN  GOTO  2880 
287G  Y^-AA<J)«BTIME<I>"<NS+1-JMY2  0  GOTO  2890 
2880  Y2-AA<J)4Y2 
2890  NEXT  J 

2900  CHI2-<BTEMP<I>-Y2>~2/Y24CHI2 
2910  NEXT  I 


2920 

2930 

2940 

2950 

2960 

2970 

2980- 

2990 

3000 

3010- 

3020 

3030 

3040 

3050 

3060 

3070 

3080 

3090 

3100 

3110 

3120 

3130 

3140 

3150 

3160 

3170 

3180 

3190 

3200 

3210 

3220 

3230 

3240 

3250 

3260 

3270 

3280 

3290 

3300 

3310 

3320 

3330 

3340 

3350 

3360 

3370 

3380 

3390 

3400 

3410 

3420 

3430 

3440 

3450 

3460 

3470 

3480 

3490 

3500 

3510 

3520 


DISP  "OBSERVED  CHI2-  ",CHI2,"  WITH  ",N+NS+I,"  D.F.” 

DISP  "DO  YOU  WANT  A  HARD  COPY  PLOT?  Y(ES)  OR  N(0)" 

INPUT  NYSft  IF  NUMSMO-1  THEN  GOTO  3040 
DISP  "TIME  WILL  BE  PLOTTED  FROM  ",XMIN,MTO",XMAX 
DISP  "TEMPERATURE  WILL  BE  PLOTTED  FROM" , YMIN,"TO*\YMAX 
DISP  "IF  YOU  ARE  SATISFIED  PRESS  RETURN" 

INPUT  NYLIMS 

IF  NYLIMS-""  THEN  GOTO  3040 

DISP  "ENTER  TIME  RANGE  TO  BE  PLOTTED,  THIN,TMAX" 

INPUT  XMIN.XMAX 

DISP  "ENTER  TEMPERATURE  RANGE  TO  BE  PLOTTED" 

INPUT  YMIN.YMAX 

IF  NTEST>-  1  THEN  GOTO  3340 

CLEAR 

DISP  "CALCULATING  PLOT  SCALING  PARAMETERS  "&FN I V IDS < "PLEASE  WAIT") 
AXMAX-INT  (LGT  <XMAX>> 

AYMAX-INT  (LGT  <YMAX>> 

XF-1 0'-AXMAX  $  YF-10~AYMAX  ft  DELX-.01  6'  DELY-.01 
XPWRS-VALS  < AXMAX)  ft  YPWRS-VALS  (AYMAX) 

!  XMIN-0  6  YMIN-0 

IF  <XMAX-XMIN)/XF>-  .1  THEN  DELX-.02 

IF  (XMAX-XMIN) /XF>-  .3  THEN  DELX-.05 

IF  (XMAX-XMIN) /XF>-  .5  THEN  DELX-.1 

IF  (XMAX-XMIN) /XF>-  1  THEN  DELX-.2 

IF  (XMAX-XMIN) /XF>-  3  THEN  DELX-.5 

IF  <XMAX-XMIN)/XF>-  5  THEN  DELX-1 

IF  ( YMAX-YMIN)/YF>-  .1  THEN  DELY-.02 

IF  (YMAX-YMIN)/YF>«  .3  THEN  DELY-. 05 

IF  (YMAX-YMIN)/YF>«  .5  THEN  DELY-.1 

IF  <YMAX-YMIN)/YF>-  1  THEN  DELY-.2 

IF  ( YMAX-YMIN)/YF>-  3  THEN  DELY-. 5 

IF  < YMAX-YMIN)/YF>-  5  THEN  DELY-1 

DISP  "THE  X-AXIS  TICK  SPACING  WILL  BE”,DELX 

DISP  "IF  YOU  ARE  SATISFIED  PRESS"AFNIVIDS("ENDLINE") 

INPUT  NTKSft  IF  NTKS=""  THEN  GOTO  3280 

DISP  "ENTER  X-AXIS  TICK  SPACING"  ft  INPUT  DELX 

DISP  "THE  Y-AXIS  TICK  SPACING  WILL  BE", DEL Y 

DISP  "IF  YOU  ARE  SATISFIED  PRESS"AFNIVIDS("ENDLINE"> 

INPUT  NTKSft  IF  NTKS-""  THEN  GOTO  3320 

DISP  "ENTER  Y-AXIS  TICK  SPACING"  ft  INPUT  DELY 

XNUM-CEIL  ( ( XMAX-XMIN) /( XF*DELX >  >  ■„ 

YNUM-CEIL  < ( YMAX-YMIN > /( YF*DEL Y ) ) 

IF  NY$-"N"  THEN  GOTO  3380 
PLOTTER  IS  505 

GRAPH  ft  IF  NUMPLT-0  THEN  GCLEAR  ft  DEG 
GOTO  3420 
PLOTTER  IS  1 
GRAPH 

f  IF  NUMPLT#0  AND  NTEST>0  THEN  GOTO  3110 

IF  NUMPLT-0  THEN  GCLEAR 

LOCATE  30, RATIO  •100-10,20,95 

t  IF  NUMPLT#0  AND  NTEST>0  THEN  GOTO  3110 

SCALE  XMIN/XF , XMIN/XF+XNUM-DELX , YMIN/YF , YMIN/YF* YNUM*DELY 

IF  NTEST>0  AND  NUMPLT#0  THEN  GOTO  3670 

FXD  2 

LGRID  DELX, DELY, XMIN/XF, YMIN/YF 
AXS- (XMAX-XMIN) /(XF*10> 

AYS-( YMAX-YMIN )/(YF«1 0) 

MOVE  (XMAX*XM1N)/(2«XF> , YMIN/YF-2»AYS  ft  LORG  6 
LABEL  USING  "K"  ;  "TIME  (SECS)  x  1E"AXPWR$ 

MOVE  XMIN/XF-2*AXS,( YMAX-YMIN) /(2»YF> 


3530  LGRG  4  9  DEG  &  LDIR  90 

3540  LABEL  USING  "K”  ;  "TEMP.<DEGS  C>  x  1E"AYPWR$ 

3550  LDIR  0 

3560  IF  NPROBE-2  THEN  GOTO  3630 
3570  LINE  TYPE  1 

3580  MOVE  AT IME  <  JM IN  > /XF , ATEMP  <  JMIN  > / YF 
3590  FOR  I- JMIN  TO  JMAX 
3600  DRAW  ATIME(I)/XF ,ATEMP(I>/YF 
3610  NEXT  I 

362fr  IF  NPROBE-1  THEN  GOTO  3660 
3630  LINE  TYPE  5 

3640  MOVE  BTIME(KMIN>/XF,BTEMP<KMIN)/YF 

3650  FOR  I-KMIN  TO  KMAX  $  DRAW  BTIME(I>/XF,BTEMP<I)/YF  «  NEXT  I 
3660  IF  NYSL$-"N”  OR  NYSL$-"NO"  THEN  GOTO  3720 
3670  LINE  TYPE  4 

3680  MOVE  BTIME(801 )/XF,BTEMP<801 )/YF 

3690  FOR  1-801  TO  1000 

3700  DRAW  BTIME<I)/XF ,BTEMP<I)/YF 

3710  NEXT  I 

3720  PEN  UP 

3730  MOVE  XMIN/XF, YMIN/YF 

3740  IF  SAR-0  THEN  GOTO  3780 

3750  MOVE  XMIN/XF+4 .5*AXS, YMAX/YF+ . 1*AYS 

3760  LORG  4 

3770  LABEL  "SAR",SAR,"W/I<g" 

3780  NTEST-NTEST+1 

3785  DISP  “PRESS  "AFNIVIDS<"CONTINUE">&"TQ  CALCULATE  SAR“ 

3790  PAUSE 
3800  KEY  LABEL 
3810  RETURN 

3820  Find:  !  ***  SUBROUTINE  TO  FIND  POINT  WHERE  HEATING  STARTS 
3830  DISP  “ENTER  A  TEST  RANGE  J1 , J2“  6  INPUT  J1,J2 
3840  DISP  "TIME" ."TEMP", “DT/Dt“ 

3850  FOR  I-J1  TO  J2-1 

3860  DISP  ATIME< I) , ATEMP< I) , (ATEMP < 1+1 )-ATEMP<I> >/<ATIME(I+1 )-ATIME(I) ) 
3870  NEXT  I 

3880  DISP  "TRY  ANOTHER  RANGE  -  Y<ES)  OR  N<0>"  6  INPUT  NYR$ 

3890  IF  NYR$="Y"  OR  NYR£-"Y ES"  THEN  GOTO  3830 
3900  DISP  "ENTER  JTURN"  8  INPUT  JTURN8  RETURN 
3910  CHAIN  "DUMP" 

3920  END 


PROGRAM  FOR  DATA  ANALYSIS,  CURVE  SMOOTHING  AND  GRAPHICS  -  version  B 


to  !  . . * . * . ** . 

20  !  * 

30  !  *  PROGRAM  TO  READ  TEMPERATURE  READINGS  TAKEN  FROM  * 

AO  !  *  TWO  NARDA  PROBES  AND  STORED  ON  DISC.  PLOT  DATA  * 

50  !.  *  IN  VARIOUS  FORMS  AND  DO  DOSIMETRY  CURVE  FITTING  * 

60  !  *  * 

70  !  *  WRITTEN  BY  CHRISTOPHER  C  DAVIS.  DECEMBER.  1986  * 

80  !  *  * 

90  ♦'  **************** . ****** . *********** . 

100  OPTION  BASE  1 

HO  MASS  STORAGE  IS  ":D500" 

120  DIM  ATEMPl  1 00 0 > . ATIMEt 1 000) .BTEMPl 1000) ,BTIME< iOOO) ,PQNT<50> ,P0FFT(50> 
130  DIM  ANUM(20> ,BNUM(20> ,IPDN<20> ,IP0FF<20> .DUMMYS! 501 
140  !  ***  FUNCTION  TO  GIVE  INVERSE  VIDEO  DISPLAYS 
150  DEF  FNIVIDS(DUMMYS) 

160  LENGTH-LEN  (DUMMYS) 

170  FOR  II-1  TO  LENGTH 

180  DUMMYS t II , II ] -CHRS  (NUM  (DUMMYS! II , III >+t28> 

190  FNIVIDS- DUMMYS 
200  NEXT  II 
210  FN  END 
220  CLEAR 

230  DISP  FNIVIDSC-DATA  DATE?**) 

240  DISP  FNIVIDSC’DATE  FORMAT" > 

250  DISP  FNIVIDSC'DAY  MONTH , YEAR") 

260  INPUT  EXPDATES 

270  DISP  FNIVIDSC’ENTER  SAMPLE  ID") 

280  INPUT  SAMPLES 

290  EXPDATAS-EXPDATES&SAMPLES 

300  ASSIGN*  1  TO  EXPDATAS&" :D501 " 

310  NUMSMO-O  v  NTEST-0  9  SAR-0 

320  IATOT-O  $  IBTOT-O  9  IONT-O  $  IOFFT-O  S  IAt-1  d  IBl-l  »  ION 1  - 1  *  IOFF1-I 

330  READ*  1  ;  IATOT  !  READ  NO  OF  PROBE  A  READINGS 

340  FOR  I-IA1  TO  IATOT  $  READ*  1  :  ATEMP(I) .ATIME<I)$  NEXT  I  !  PROBE  A  DATA 

350  READ*  1  ;  IBTOT  !  READ  NO  OF  PROBE  B  READINGS 

360  FOR  I-IB1  TO  IBTOT  9  READ*  1  ;  BTEMP< I ) ,BTIME< I)$  NEXT  I  !  PROBE  B  DATA 

370  READ*  1  ;  10NT  !  NO  OF  TIMES  MW  PWR  TURNED  ON 

380  FOR  I-IONI  TO  IONT  9  READ*  1  ;  PQNT(I>$  NEXT  I  f  ON  TIMES 

390  READ*  1  ;  IOFFT  !  NO  OF  TIMES  MW  PWR  TURNED  OFF 

4C0  FOR  I-IOFF1  TO  IOFFT  *  READ*  1  ;  POFFT(I>$  NEXT  I  !  OFF  TIMES 

410  IA1-IA1+IATOT  $  IB1 -IB1 +IBT0T  9  ION1-ION1+IONT  6  IOFF1 -IDFFl +IOFFT 

420  ASSIGN*  1  TO  * 

430  !  PURGE  EXPDATAS&" :D501" 

440  DISP  "TOTAL  PROBE  A  READINGS  -  ", IATOT 

450  !  FOR  1-1  TO  IATOT  *  DISP  ATEMP< I) , ATIME(I) , I  9  NEXT  I 

460  DISP  "TOTAL  PROBE  B  READINGS  -  ".I8T0T 

470  !  FOR  1-1  TO  IBTOT  6  DISP  BTEMP(I)  , BT IME < I >  ,1  6  NEXT  I 

4601  DISP  "MW  POWER  TURNED  ON  *\IONT,"  TIMES" 

490  DISP  "MW  POWER  TURNED  OFF  ".IOFFT,"  TIMES" 

500  DISP  "THE  MW  PWR  WAS  TURNED  ON  AT  THE  FOLLOWING  TIMES" 

510  FOR  1-1  TO  IONT  %  DISP  PONT(I)  «  NEXT  I 
520  DISP  "Sample  ID  is  " ,FNIVID$(SAMPLE$> 

530  DISP  "Date  of  exposure  was  " .FNIVIDS(EXPDATES) 

540  1 

550  ! 

560  ! 

570  1 

580  ! 

590  ! 

600  ! 


SECTION  FOR  SMOOTHING 
AND  PLOTTING  OF  EXPERIMENTAL  DATA 

WRITTEN  BY  Christopher  C.  Davis,  June  1984, 
MODIFIED  DECEMBER  1986 


620 

630 

£40 

G 

650 

660 


DIM  AA<5),SN(21,21),R<20) ,HEAT <5> ,CQ0L(5>  !  ARRAYS  USED  IN  CURVE  SMOOTHIN 


NYSLS-"N" 

IF  ATIMEl3AT0T>>BTIME(IBT0T> 


>  !  SET  X  AXIS  MAX  TO  MAX  OBSERVATION  TIME 


THEN  XMAX-ATIME(IATQT)  ELSE  XMAX-BTIMEdBTOT 


670 
680 
690 
700 
710 
720 
730 
740 
750 
760 
770 
780 
790 
800 
810 
OR  3") 
820 
830 
840 
850 
860 
870 
880 
890 
900 
910 
920 
930 
940 
950 
960 
970 
980 
990 
1000 
1010 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1110 
1120 
1130 
1140 
1150 
1160 
1170 
1180 
1190 


XMIN-0  !  INITIALLY  SET  X  ORIGIN  TO  0 
JMIN-1  6  KMIN-1 

DISP  "TEMPERATURE  MEASUREMENT  INTERVAL  WAS",XMAX 
YMIN-10‘20  6  YMAX—  <10~20> 

FIND  MAXIMA  AND  MINIMA  OF  DATA 
FOR  1-1  TO  IATOT 

IF  ATEMP< I)>YMAX  THEN  YMAX-ATEMPd) 

IF  ATEMP<IXYMIN  THEN  YMIN«ATEMP<I> 

NEXT  I 

FOR  1-1  TO  IBTOT 

IF  BTEMP< I)>YMAX  THEN  YMAX-BTEMP< I> 

IF  BTEMP<IXYMIN  THEN  YMIN-BTEMPd) 

NEXT  I 

NUMPLT-0  g  IF  NTEST#0  THEN  CLEAR 

DISP  "DO  YOU  WISH  TO  PLOT  PROBE  A,  PROBE  B,  OR  BOTH? 


*&FNrviDS<"ENTER  1,2 


DISP  "YOU  CAN  ONLY  CALCULATE  SAR  USING  PROBE  A" 

INPUT  NPROBE 

IF  NPROBE- 1  THEN  GOTO  860 
GOTO  870 

FOR  I-J  TO  IATOT  S  BTIMEdXATIMEd)  *  BTEMP ( I > - ATEMP Cl >  *  NEXT  I 

IF  NTEST-0  THEN  GOTO  890 

NYSLS-"Y" 

DISP  "ENTER  TIME  REGION  YOU  WISH  TO  ANALYZE" 

INPUT  TMIN.TMAX 

FOR  1-1  TO  IATOT  f  FIND  RANGE  OF  TEMP  ARRAY  TO  BE  PLOTTED 
IF  TMIN>ATIMEd>  THEN  JMIN-I 
NEXT  I 

FOR  1-1  TO  IBTOT  8  IF  TMIN>BTIME<I>  THEN  KMIN-I 
NEXT  I 

IF  ATIME<  JMINXBTIME<KMIN>  THEN  TMIN-ATIME1 JMIN)  ELSE  TMIN-BTIME<KMIN> 
FOR  1-1  TO  IATOT  #  IF  TMAX>ATIME(I)  THEN  JMAX-I 
NEXT  I 

FOR  I- 1  TO  IBTOT  *  IF  TMAX>BTIMEU>  THEN  KMAX-I 
NEXT  I 

IF  ATIME< JMAX)>BTIME<KNAX>  THEN  TMAX-ATIME( JMAX)  ELSE  TMAX-BTIMECKMAX) 
IF  NTEST-0  THEN  GOTO  1250 
IF  NYSL$-"N"  OR  NYSL$-"NO"  THEN  GOTO  1250 
!  ***  SECTION  FOR  SELECTION  OF  CURVE  REGION  FOR  SMOOTHING 
TEMAX--1.E20  #  TEMIN-1.E20 

FOR  I- JMIN  TO  JMAX  #  IF  ATEMPdXTEMIN  THEN  TEMIN-ATEMPd) 

IF  flTEhP(I)-TEMIN  THEN  JTURN-I  !  Find  inflection  point  on  temp  curve 

DISP  "NUMBER  OF  POINTS  ON  COOLING  SLOPE-  ", JTURN-JMIN+1 
DISP  "NUMBER  OF  POINTS  ON  HEATING  SLOPE-  ", JMAX-JTURN+1 
DISP  "IS  THIS  SATISFACTORY?  Y<ES>  OR  N(0>"  »  INPUT  NYSATS 
IF  NYSAT$-"N"  OR  NYSAT$-"NO"  THEN  GOTO  890 
IF  NYSLS-"N"  OR  NYSL$«"NO"  THEN  GOTO  1250 
CLEAR  ft  DISP  "COOLING  REGION  POINTS", JMIN, JTURN 
DISP  "HEATING  REGION  POINTS", JTURN, JMAX 
ON  KEY#  1, "COOLING"  GOTO  1210 
ON  KEY#  2, "HEATING"  GOTO  1230 
DISP  FNIVID$<"SELECT  OPTION")  6  DISP  •  DISP 
„ _ KEY  LABEL 


200  GOTO  1200 

210  NSM-l  9  DISP 

r  JMINt JTURN 
(220  GOTO  1240 

1230  NSM-2  9  DISP 

fr  JTURN, JMAX 


•ENTER  POINTS  RANGE  FOR  COOLING  SLOPE"  C«  DISP  9  DISP  S>  INPU 
•ENTER  POINTS  RANGE  FOR  HEATING  SLOPE"  9  DISP  9  DISP  9  INPU 


1240 

1250 

1260 

1270 

1280 

1290 

300 

1310 

1320 

1330 

1340 


NUMPLT-1  9  GOSUB  Smooth 
GOSUB  Aplot 

ON  KEY*  1, "AGAIN"  GOTO  1380 
ON  KEY*  2, "DOSE"  GOTO  1430 
ON  KEY*  3, "PLOT"  GOTO  1440 
ON  KEY*  4 ."SCRDMP"  GOTO  3530 
ON  KEY*  5, "EXIT"  GOTO  3540 
CLEAR  9  KEY  LABEL 

DISP  "SELECT  OPTION  -  "aFNIVID$<"DOSE">a"  ALLOWS  CALCULATION  OF  SAR" 
DISP  "  -  "&FNIVIDS( "AGAIN" >&"  RE-READS  DATA  FROM  DISC" 

DISP  "  -  "aFNIVIDS<"PLQT">a"  RE-PLOTS  LAST  PLOT  -  USE  WITH 


1350  DISP  "  -  "aFNIVIDS("SCRDMP">a"  DUMPS  GRAPHICS  TO  PRINTER" 

1360  DISP  "  -  "aFNIVIDS("EXIT">a"  TERMINATES  EXECUTION" 

1370  GOTO  1370 

1380  CLEAR 

1390  DISP  "IF  YOU  CONTINUE  DATA  WILL  BE  RE-READ  FROM  DISC" 

1400  DISP  "PRESS  "4FNIVIDS< "AGAIN" >a"KEY  IF  YOU  WISH  TO  DO  THIS" 

1410  ON  KEY*  1, "AGAIN"  GOTO  300  $  KEY  LABEL 
1420  GOTO  1420 
1430  GOTO  860 

1440  IF  NUMSMO-O  THEN  GOTO  1470 

1450  GOSUB  Aplot  9  GOTO  1370 

1460  GOTO  1260 

1470  DISP  "OPTION  NOT  AVAILABLE"  9  GOTO  300 

1480  Smooth: 

H30  ♦  . . . 

1500  1  * 

1510  !  -  CURVE  SMOOTHING  SECTION 

1520  !  * 

1530  ! . . . . 

1540  DISP  "CURVE  FITTING  IN  PROGRESS  "aFNIVIDS ("PLEASE  WAIT") 

1550  NUMSMO-1 

1560  IF  NSM-2  THEN  GOTO  1600 

1570  N-JTURN-JMIN+1 

1580  FOR  I-t  TO  N  9  BTIME<I)-BTIME< JMIN+I-1 >  9  BTEMPU >-BTEMP< JHIN+I-1 >  %  NEX 
T  I 

1590  GOTO  1630 

1600  N- JMAX- JTURN+1 

1610  FOR  1-1  TO  JMAX- JTURN+1  9  BTIMEC I>-BTIME( JTURN+I-1 )  9  BTEMP<I>-BTEMP( JTU 

RN+-I-1) 


1620  NEXT  I 

1630  FOR  1-1  TO  5 

1640  R(I)-0  9  AA< I )-0  9  FOR  J-1  TO  21  9  SN<r,J>-0  9  NEXT  J  9  NEXT  I 

1650  DISP  "WHAT  ORDER  OF  POLYNOMIAL  DO  YOU  WISH  TO  FIT  -  1 -LINEAR. . ." 

1660  INPUT  NS 

1670  FOR  J-1  TO  NS+1 

1680  FOR  I-J  TO  N 

1690  IF  J-NS+1  AND  BTlME<I)-0  THEN  GOTO  1710 

1700  R<J)-BTEMP<I>-BTIME<ir<N$+1-J)*R<J>  9  GOTO  1720 

1710  R(J>-BTEMP(I>+R(J) 

1720  NEXT  I 

1730  NEXT  J 

1740  FOR  K-1  TO  NS+1  9  FOR  J-1  TO  NS+1  *  FOR  1-1  TO  N 

1750  IF  2-NS-J-K+2-G  AND  BTIME<I>-0  THEN  GOTO  1  l/U 

1760  SN<J,K)-BTIME(I>  <2«»NS-J-K+2MSN<  J,K>  9  GOTO  1760 


780 
n  790 
1800 
1810 
1820 
1830 
1840 

[1850 
1860 
DR  J-1 
1870 
1880 
1890 
1900 
1910 
1920 
1930 
1940 
1950 
1960 
1970 
1980 
1990 
2000 
2010 
2020 
2030 
2U40 
2050 
2Ub0 
2070 
2080 
2090 
2100 
2110 
2120 
2130 
2140 

21SU 

2160 

2170 

2180 

2190 

22QQ 

2210 

2220 

2230 

2240 

2250 

2260 

2270 

2280 

L 

2290 

2300 

2310 

2220 

2330 

2340 


I'ii'iv  j  ,k/ 

NEXT  I 
NEXT  J 
NEXT  < 

GOSUB  Irwt 

FOR  1-1  TO  NS+ 1  «  FOR  J-1  TO  NSH 
AA<NS+2-I)-SN<NS+2-I , J>*R< J)+AA<NS*2-I>  9 
IF  NSM-1  THEN  DX-ATIME( JTURN>-ATIME< JMIN) 


NEXT  J  €>  NEXT  I 

ELSE  DX-AT I ME  < JMAX ) - ATIME l JTUR 


IF  NSM-1  THEN  XSTART-ATIHE< JMIN)  ELSE  XSTART-ATIME< JTURN) 

FOR  K-1  TO  200  9  I-K+800  9  BT IHE < I > - < K- 1 >*DX/20G+XSTART  6  BTLMP( 
TO  NS+1 

IF  NS+1-J-0  AND  BTIM£(I)-0  THEN  GOTO  1900 
BTEMP<I)-AA(J)*BTIHE(1)'  <NS+1-J)+8»lHPU> 

GOTO  1910 

BTEMP<I)-AA< J)+BTEMP(I> 

NEXT  J 
NEXT  K 
RETURN 
Invt : 


P-0  9 


F 


MATRIX  INVERSION  BY  GAUSS-JORDAN  ELIMINATION  WITH  FULL  PIVOTING 
SQUARE  MATRIX  OF  ORDER  NP 


DIM  KP<5) , JP<5) 

NP-NS-M 

NCOL-NP 

♦ 

!  INITIALIZE  PERMUTATION  VECTORS 

♦ 

‘  FOR  1-1  TO  NP  £  KP< I)-I  $  JP(I)-I  9  NEXT  I 
* 

!  INVERT  MATRIX 

f 

*  FOR  IR= 1  TO  NP 
f  FIND  THE  PIVOT  ELEMENT 
PIVOT-O 

FOR  I-IR  TO  NP  *  FOR  J-IR  TO  NP  9  K-KP<I>  9  L-JP(J) 
lESl-Abb  <SN<K,L>>-HBb  tPlVUi; 

IF  TEST<-  0  THEN  GOTO  2220 

IF  TEST>G  THEN  GOTO  2190 

GOTO  2220 

IK-I 

JK-J 

PIVOT-SN(K ,L) 

NEXT  J  v  NEXT  I 

! 

!  THE  PIVOT  ELEMENT  IS  SN<K,L> 

! 

1  UPDATE  PIVOT  VECTORS 
» 

K-KPUK)  i>  KP(IK)-KPdR)  9  KP(IR)-K  9  L-JP(JX)  9  JP< JK>-JP(IR>  9  JP(IR>- 


! 

!  COMPUTE  NEW  ELEMENTS  OF  PIVOTAL  ROW 
* 

j  FOR  J-1  TO  NCOL  9  SN<K , J1-SN1K , J1/PIVOT  9  NEXT  J 
!  COMPUTE  REMAINING  ELEMENTS  OF  IR'TH  STEP 


40JU 
AstOM 

2370 
2380 
2390 
2400 
2410 
2420 
2430 
2440 
2450 
2460 
2470 
*480 
2490 
2500 
2510 
2520 
2530 
2540 
2550 
2560 
2570  NLT-1 

2580  CLEAR  £  IF  NUMPLT-0  THEN  GOTO  2780 
2590  ALPHA  £  DISP  "PARAMETERS  OF  FIT"  £  DISP 
2600  FOR  1-1  TO  NP  £  DISP  AA < I >  6  NEXT  I 
2blO  IF  NSM-2  THEN  GOTO  2640 
2620  FOR  1-1  TO  NP  £  COOL<I)-AA(I)  £  NEXT  I 
*b30  GOTO  2650 

2640  FOR  1-1  TO  NP  £  HEAT(I)-AA<I>  £  NEXT  I 
2650  CHI2-0 

2660  IF  NTEST<2  THEN  GOTO  2700 
2670  SAR-4 . 2* < HEAT ( 1 ) -COOL ( 1 > ) 

2680  IF  NS-2  THEN  SAR-2-SAR-ATIME< JTURN)*4.2*(HEAT<2>-C0QL<2) > 

2690  DISP  FNIVIDS<"SAR  IS") ,SAR , "Watts  per  gTaw" 

2700  FOR  1-1  TO  N  £  Y2-0  £  FOR  J-1  TO  NS+1 
2710  IF  NS+1-J-0  AND  BTIME<I)-0  THEN  GOTO  2730 
2720  Y2-AAC J)*BTIME(I)~(NS+1-J)+Y2  £  GOTO  2740 
2/30  Y2-AA(J)+Y2 
2740  NEXT  J 

2750  CHI2-<BTEHP(I>- Y2)“2/Y2+CHI2 
2760  NEXT  I 

2770  DISP  "OBSERVED  CHI2-  ".CHI2,"  WITH  '\N*NS+1 D.F." 

2/8U  DISP  "DO  YOU  WANT  A  HARD  COPY  PL017  Y(£S)  OK  N(U)' 

2790  INPUT  NYS£  IF  NUMSMQ-1  THEN  GOTO  2890 

2800  DISP  "TIME  WILL  BE  PLOTTED  FROM  ",XMIN,"TO",XMAX 

2810  DISP  "TEMPERATURE  WILL  BE  PLOTTED  FROM", YMIN,"TO" ,YMAX 

2820  DISP  "IF  YOU  ARE  SATISFIED  PRESS  RE1UKN" 

2830  INPUT  NYLIMS 

2840  IF  NYLIMS-""  THEN  GOTO  2890 

2850  DISP  "ENTER  TIME  RANGE  TO  BE  PLOTTED,  THINJMAX" 

*8 bU  ANPUl  XMAN ,XMAX 

2870  DISP  "ENTER  TEMPERATURE  RANGE  TO  BE  PLOTTED" 

2880  INPUT  YMIN, YMHX 

2890  IF  NTEST>-  1  THEN  GOTO  3050 

2900  CLEAR 

2910  DISP  "CALCULATING  PLOT  SCALING  PARAMETERS  "aFN IV IDS ("PLEASE  WAIT") 

*8*0  axmax-in;  <lgt  <xmax)) 

2930  AYMAX-INT  (LGT  <YMAX>) 

2940  XF-10*  AXMAX  i  YF-10AYMAX 

2950  XPWRS-VAL5  <AXMAX>  £  YPWRS-VALS  <AYMAX) 

2960  !  XMIN-0  £  YMIN-O 


run  x* »  iu  nr 

IF  I-K-0  THEN  GOTO  2440 
AIL-SN< I ,L> 

SN<I,L>— (AIL/PIVOT) 

FOR  J-1  TO  NCOL 
IF  J-L-0  THEN  GOTO  2430 
SN<I , J)*SN<I , J)-AIL*SN(K , J> 
NEXT  J 
NEXT  I 

* 

! 

!  INVERT  PIVOT  ELEMENT 
! 

SN<K,L)- 1 /PIVOT  £  NEXT  IR 
RETURN 
Aplot : 


m  • 

*  PLOTTING  SUBROUTINE  * 


2970  IF  XMAX/XF>«  1  THEN  DELX-.2 
2980  IF  XMAX/XF>-  3  THEN  DELX-.5 
2990  IF  XHAX/XF>-  5  THEN  DELX-1 
3000  IF  YMAX/YF>-  1  THEN  DELY-.2 
3010  IF  YMAX/YF>-  3  THEN  DELY-.5 
3020  IF  YMAX/YF>-  5  THEN  DELY-1 
3030  XNUM-CEIL  < <XMAX-XMIN>/<XF*DELX) > 

3040  YNUM-CEIL  <<YMAX-YMIN)/<YF*D£LY>> 

3050  IF  NYS-’W  THEN  GOTO  3090 
3060  PLOTTER  IS  505 

3070  GRAPH  $  IF  NUMPLT-0  THEN  GCLEAR  @  DEG 
3U80  GOTO  3130 
3u90  PLOTTER  IS  l 
3100  GRAPH 

3110  !  IF  NUMPLT#0  AND  NTEST>0  THEN  GOTO  3110 

3120  IF  NUMPLT-0  THEN  GCLEAR 

3130  LOCATE  30, RATIO  *100-10,20,95 

3140  !  IF  NUMPLT#0  AND  NTEST>0  THEN  GOTO  3110 

3150  SCALE  XMIN/XF ,XMIN/XF+XNUM*DELX , YMIN/YF , YMIN/YF+YNUM*DELY 

3lbU  if  NILST>0  ANU  NUMPL I *U  THEN  GUiU  33BU 

3170  FXD  1 

3180  LGRID  DELX.DELY, XMIN/XF, YMIN/TF 
3190  AXS«(XMAX-XMIN)/(XF*10) 

3200  AY$-( YMAX-YMIN)/(YF*1 0) 

3210  MOVE  XMIN/XF+4 .5*AXS , YMIN/YF-2*AYS  3  LORG  6 
3220  LABEL  USING  "K"  ;  "TIME  <SECS>  x  1EM*XPWR$ 

^23U  MUVt  XMIN/XF -2* AXS , YMIN/ YF+4 .5*AtS 
3240  LORG  4  3  DEG  3  LDIR  90 

3250  LABEL  USING  "K"  ;  "TEMP . (DEGS  C)  x  1E"&YPWR$ 

3260  LDIR  0 

3270  IF  NPR0BE-2  THEN  GOTO  3340 
3280  LINE  TYPE  1 

I  3290  MOVE  ATII1E<  JNIN)/XF,  ATEMP<  JMINl/YF 
33UO  FOR  I-JMIN  TO  JMAX 
3310  DRAW  ATIMEi I)/XF , ATEMP(I>/YF 
3320  NEXT  I 

3330  IF  NPROBE-!  THEN  GOTO  3370 
3340  LINE  TYPE  5 

3350  MOVE  BTIME<KMIN)/XF,BTEMP(KMIN)/YF 

3360  FOR  I-KMIN  TO  KMAX  6  DRAW  BTIME(I)/XF,BTEMP(I)/YF  %  NEXT  I 
3370  IF  NYSL$-"N"  OR  NYSL$-"NO"  THEN  GOTO  3430 
338U  LINE  TYPE  4 

3390  MOVE  BTIME<801 )/XF,BTEMP<801 >/YF 

3400  FOR  1-801  TO  1000 

3410  DRAW  BTIME( I)/XF ,BTEMP< I )/YF 

3420  NEXT  I 

3430  PEN  UP 

3440  MOVE  XMIN/XF, YMIN/YF 

3450  IF  SAR-0  THEN  GOTO  3490 

3460  MOVE  XMIN/XF+4.5*AXS, YMAX/YF*. 1*AYS 

3470  LORG  4 

3480  LABEL  "SAR" ,SAR , "Watts  per  gran" 

3490  NTEST-NTEST+1 
3500  PAUSE 
3510  KEY  LABEL 
3520  RETURN 
3530  CHAIN  "DUMP" 

35*40  END 
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INTRODUCTION 

Dielectric  heating  is  a  recognized  mechanism  for  the  induction  of  radiofrequency  (RF) 
bioeffects  considered  to  be  thermal  in  origin  •  that  b  attributable  to  temperature  increases. 
However,  apart  from  RF-beat  cell  killing  (e.g.  Sapareto  et  al.,1982,  Chang  et  al.,1987)  very 
little  information  b  available  on  quantitative  relationships  between  RF  bioeffects  at  the 
cellular  and  molecular  level  and  temperature  profiles  over  time  or  thermal  dosage  (TD).  The 
determination  of  RF  heating/cooling  curves  can  be  used  to  compute  the  specific  absorption 
rate  (SAR)  and  the  specific  absorption  (SA)  in  an  exposed  sample  (Stuchly  and  Stuchly, 
1986).  Biological  variables  and  exposure  conditions  can  be  controlled  in  experiments  with  in 
vitro  systems  to  a  degree  not  achievable  in  vivo.  An  attractive  model  to  study  RF  bioeffects 
b  the  transformation  of  lymphocytes  in  vitro,  provided  the  biological  variables  can  be  related  ; 
to  dosimetric  quantities  that  characterize  exposure  (Czerski,1975;  Budd  and  Czerski,1985).  j 
To  accomplish  thb,  we  designed  an  exposure  system  with  provbion  for  real-time  temperature  : 
monitoring  with  RF-field  non-perturbing  temperature  probes.  The  exposure  system,  which 
will  be  described  in  detail  later,  has  multiple  sample  chambers.  To  allow  on-line  thermometry 
and  dosimetry,  one  of  these  chambers  b  used  as  a  site  for  a  non-perturbing  temperature  probe. 
From  the  temperature  (T) /time(t)  hbtory  of  the  sample  chamber,  an  exposure  dosage  can  be 
determined  and  described  in  terms  of  TD,  SAR,  SA,  or  electric  field  strength.  The  system  can , 
be  used  for  studies  of  RF-bioeffects  in  any  tissue  culture  cell  line  or  other  in  vitro  biological  | 
sample.  The  monitoring  component  can  be  applied  to  the  study  of  ‘emperature-dependent  i 
effects,  irrespective  of  the  modality  used  for  heating.  | 

EXPOSURE  SYSTEM  1 


The  output  from  a  CW  Hewlett-Packard  Model  8616A  oscillator  feeds  a  Hughes  Model  1177H 
TWT  amplifier.  Pulsed  or  amplitude-modulated  waveforms  are  obtained  through  the  use  of 
a  Hewlett-Packard  Model  8403A  p-i-n  modulator  and  driver.  For  high  peak  power,  but  low 
duty-cycle,  pulsed  exposures  the  signal  can  be  further  amplified  with  a  lkW  amplifier  (MCL, 
La  Grange,  Ill.)  The  amplified  signal  enters  a  shorted  section  of  S-band  rectangular  waveguide 
300  mm  long  through  a  matched  coaxial  feed.  The  waveguide  sample  holder  receives  its 
power  through  an  isolator  (not  shown  specifically  in  Fig.  1)  and  a  dual-directional  coupler 
that  allows  forward  and  reflected  powers  to  be  measured  with  Hewlett-Packard  Model  432A 
power  meters,  and  waveforms  to  be  monitored.  A  sample  holder  for  biological  specimens 
is  supported  centrally  in  the  waveguide  in  a  block  of  low-density  polystyrene  that  fills  the 
cross-section  of  the  guide.  An  identical  waveguide  section  with  an  identical  plastic  sample 
holder  is  used  for  sham  exposures.  During  microwave  exposures,  both  waveguide  assemblies 
were  housed  in  a  COj  tissue  culture  incubator  thermostated  ai  37.0°C.  Various  sample 
holders  were  tested:  a  4-chamber  plastic  tissue  culture  slide  (Miles  Scientific  #  4804)  was 
found  satisfactory  in  our  application.  Each  chamber  holds  1  ml.  A  miniature  thermistor 
temperature  probe  enters  each  waveguide  section  through  a  small  hole  in  the  top  shorting 
plate.  The  hole  is  drilled  near  the  guide  wall,  at  the  center  of  the  shorter  dimension  -  this  is 
a  low  electric  field  point  and  the  hole  produces  minimal  disturbance.  The  active  end  of  the 
thermistor  probe  enters  the  sample  in  one  of  the  4  chambers  of  the  tissue  culture  dish.  We 
could  move  the  probe  from  chamber  to  chamber  to  check  exposure  uniformity.  Two  types  of 
non-field  perturbing  probe  have  been  used  in  this  way:  A  Narda  Model  801  IB  non-perturbing 
double  temperature  probe  or  two  Vitek  Model  101  probes.  The  temperatures  in  both  the 
exposed  and  sham-exposed  samples  are  recorded  continuously  during  an  experimental  run. 
The  temperature  probes  are  connected  to  a  Hewlett-Packard  relay  activator,  Model  59306A, 
which  is  itself  connected  to  a  Keithley  Model  192  DVM.  Both  the  relay  activator  and  DVM 
are  under  the  control  of  a  Hewlett-Packard  Model  86  desk-top  computer  through  the  HPEB 
(IEEE-488  bus)(Fig.2). 
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Fig. 2.  Schematic  diagram  of  multiple-chamber  sample  holder  with  non- 
perturbing  temperature  measurement  and  automated  dosimetry  system. 


Under  computer  control,  the  temperatures  of  exposed  and  sham-exposed  samples  are 
recorded  sequentially  and  stored  in  memory.  Since  the  temperature  of  sham-exposed  samples  , 
is  essentially  constant,  it  is  monitored  less  frequently  than  the  temperature  of  the  exposed 
samples.  A  typical  experimental  protocol  would  involve  10  temperature  readings  of  the  ex-  i 
posed  sample  for  every  one  reading  of  the  sham-exposed.  For  conventional  thermal  exposures, 
the  temperature  can  be  elevated  by  increasing  the  temperature  of  the  incubator. 

TEMPERATURE/TIME  ANALYSIS 

During  a  typical  experiment,  the  temperature  is  recorded  at  regular  intervals,  with  a' 
minimum  measurement  interval  <  Is,  before  the  beginning  of,  during,  and  after  the  exposure. < 
The  *on”  and  “off”  times  of  exposure  are  recorded  on  the  computer  by  the  operator  using  a  j 
“soft- key”  interrupt  capability.  At  the  conclusion  of  a  run,  the  T(t)  behavior  is  analysed  to  I 
determine  SAR.  I 

i 

If  the  SAR  is  S  (Wkg-1)  and  the  specific  heat  of  the  sample  is  C(Jkg**1K~1),  then  the 
SAK  can  be  determined  from  the  change  in  heating  rate  when  microwave  power  is  applied. 


In  practice,  this  analysis  is  performed  automatically.  At  the  end  of  a  run,  the  point  of 
inflection,  or  turning  point,  in  T(t)  is  found  numerically.  The  temperature  profile  to  the  left 
and  right  of  this  point  is  fitted  by  least-squares  to  a  linear  or  quadratic  function  and  the 
change  in  slope  at  the  beginning  of  the  exposure  period  yields  the  SAR.  Figures  3  and  4  are 
examples  of  such  a  procedure. 


Fig-3.  Temperature/time  history  of  an 
exposed  sample  showing  linear  fits  to  cool¬ 
ing  and  heating  portions  of  curve  for  dosime 
try. 


Fig.4.  Temperature/time  history  of  ex¬ 
posed  (upper  curve)  and  sham-exposed 
samples  (iower  curve).  The  dosimetry 
for  the  exposed  sample  has  been  deter¬ 
mined  from  a  quadratic  fit  to  the  heating 
curve  during  microwave  exposure  and  a 
linear  fit  to  the  equilibrated  portion  of 
the  curve  prior  to  exposure. 


Figure  3  shows  the  temperature / time  history  of  a  sample  that  was  cooling  prior  to  the  start 
of  microwave  exposure.  The  linear  fits  to  the  cooling  and  heating  portions  of  the  curve  yield 
the  SAR.  Figure  4.  shows  the  heating  of  a  sampled  that  was  quite  well  equilibrated  before 
exposure.  Because  the  signal/noise  ratio  is  high  the  SAR  was  determined  by  a  quadratic  j 
fit  to  the  heating  portion  of  the  curve.  The  quality  of  the  fit  near  the  turning  point  can  j 
always  be  examined  to  determine  whether  the  slope  value  is  realistic.  However, unless  the  i 
temperature/time  profile  has  very  low  noise,  it  is  generally  better  to  use  the  linear  fit  to  T(t) 
near  the  turning  point  or  a  biased  estimate  of  the  slope  can  result.  Figure  4  also  shows  the 
temperature/time  history  of  the  sham-exposed  sample  to  illustrate  its  temperature  stability. 

RESULTS  AND  DISCUSSION  “ 

The  advantages  of  the  procedure  described  above  are  severalfold.  The  sample  need  not  : 
be  equilibrated  before  SAR  is  determined.  Repeated  determination  of  SAR  can  easily  and  . 
quickly  be  made,  which  allows  the  SAR  uniformity  from  one  sample  chamber  to  another  : 
to  be  determined.  Exact  knowledge  of  the  microwave  power  is  not  required:  reproducible 
exposures  can  be  made  at  known  SAR  by  using  the  measured  forward  power  corresponding 
to  a  given  SAR  measurement,  provided  the  experimental  arrangement  is  not  altered  between  , 
exposures.  Dosimetry  is  not  affected  by  other  losses  in  the  system.  The  waveguide  exposure  1 
system  with  four  sample  chambers  was  found  to  give  an  SAR  uniformity  from  one  chamber 
to  the  other  within  10%.  Thus,  experimental  samples  can  be  exposed  and  examined  in  : 
triplicate,  the  fourth  chamber  being  used  for  the  insertion  of  the  temperature  probe.  I 

The  system  is  biocompatible  and,  depending  on  cell  line,  cell  density,  and  medium,  allows 
continuous  exposures  of  several  days  duration.  The  biocompatibility  of  this  arrangement  was 
proven  by  studying  the  growth  of  human  lymphocytes  under  various  conditions  within  the 
waveguide.  As  well  as  providing  biocompatible  exposure  conditions,  the  system  is  flexible: 
it  allows  the  exposure  of  tissue  culture  cells  growing  in  suspension  or  in  monolayers,  and  can 
easily  be  used  with  cell-free  samples.  The  exposure  system  is  relatively  simple,  and  can  be 
assembled  from  off-the-shelf  components.  Its  capabilities  can  be  expanded  by  introducing  . 
computer  control  of  RF  power  input  based  on  feedback  from  temperature  measurements. 
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ABSTRACT 


There  nt  contradictory  reports  in  the  literature  on  the  effects  of  micro¬ 
wave  exposure  on  transformation  of  lyaphocytes.  To  resolve  this  question, 
human  lymphocytes  have  been  exposed  in  vitro  to  2450  MH*  radiation  under 
carefully  controlled  conditions  at  various  power  levels  and  for  periods  up 
to  120  hours.  Lymphocyte  tranformation  in  exposed  samples  was  compared  to 
control  samples,  sham-exposed  or  incubated  in  a  conventional  tissue  culture 
Incubator.  Heparinized  peripheral  blood  samples  were  obtained  by  venipunc¬ 
ture  from  healthy  volunteers,  lymphocytes  were  separated  by  gradient  cen¬ 
trifugation  (Lymphocyte  Separation  Medium,  Litton  Bionetics)  and  Incubated 
in  chromosome  medium  1A  without  phytohemagglutinin  (G 1  Deo)  at  a  standard 
concentration  of  10  cells/ml  medium.  Using  the  automated  dosimetry  micro- 
wave  exposure  system  described  In  a  companion  paper  (Joyner  et  al.,  this 
Congress)  the  temperature  profile  over  time  was  recorded  for  each  sample, 
and  average  and  peak  specific  absorption  rates  (SAR),  specific  absorption 
(SA ) ,  and  thermal  dosage  (TO)  were  calculated.  At  harvest,  cells  were 
counted  and  their  viability  was  tested  using  neutral  red,  Janus  green,  and 
trypan  blue  staining.  Standardized  cytologfcal  preparation  ware  made  using 
a  "Cytospin  II*  (Shandon)  cytocentr ifuge,  air-dried,  fixed  in  methanol  and 
stained  with  a  combination  of  Wright  and  Biemsa  stains.  The  percentages  of 
untransformed  (smell)  lymphocytes,  intermediate  and  lymphoblastoid  cells 
weie  determined,  based  on  morphological  features,  quantitated  using  an 
"Optomax*  image  analyzer.  The  correlation  of  lymphocyte  transformation 
with  SAR,  SA,  and  TO  was  examined.  The  implications  of  the  results  for  a 
“thermal*  vs  *nonthermal"  mechanism  of  interaction  will  be  discussed. 
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Maryland,  Callage  run.  H)  207*2,  May*  L.  Mdoord. 
Mlatieil  Deeltk,  >«W  *m  Drue  Adaj.nl  at  rat  too , 
a.  IB  20*37. 

«•*  aystwe  Cm  aaitarini  Ut  SOSU  Chang**  in 
cm  nmi  U  active  biological  uapiM  at  a 
a*  la  baaed  an  a  aaapaar  nwwUaa  laiUtt-IMl- 
lyser,  and  uses  the  "open-probe"  technique1  with 
*  liquid  at  the  prate  tip.  Thla  frequency  doaeln 
a  iMIt  ta  26.9CSx.  la  much  aara  accurate  than 
ng  tidier  hardware2.  Ha  will  daaonatrate  a 
an  differences  in  dielectric  propertlaa  af  lata 
alea.  Thla  la  accowpliahed  by  making  no  hardware 
an  waaaurawenta .  The  ayatwa  ha*  the  additional 
blah  dearaa  of  bloeospatibillty  and  ninlnal 
•■plea .  Only  saall  quantitiea  (lOO-SOOwl)  of 
plat#  spectral  scan  ecu  be  accomplished  in  late 
awa  good  tlwe  reaelutian  in  the  monitoring  of 
tian.  *e.C.  Davit  and  C-S.  Lea,  Proc.  13th  Arm. 
.ladalpbia.  la..  March  12-13,  1*87,  pp  322-323. 
■ad  M.A.  Obeli.  Btetfva.  J.,  32.  *21-423  (1*87). 


am  inivwg  pen  modulated  so-tai  radiaticr. 

D.K.  Reoea  ,  and  w.t.  Joints.  gaalth  Li  fact  a 
m.  tel.  Ik.,  RC  27711. 

IBM.  laaoi  that  in  vitro  exposure  of  chick 
radiofrequency  (XT)  radiation,  eeplituda  nodu- 
anhanca  the  afflux  af  celciun  tana  only  within 
■to  a  perming  frew  It*  to  1.87  aM/cn2,  and  tha 
No-effect  fault*  accurrad  at  0.37,  0,7},  2.17 
ted  and  eatwHad  thaee  raaulta  by  teaLirq  nether 
iB  tha  rant*  of  the  previous  study.  Percbrain 
tfeana  wets  labeled  In  vitro  with  radioactive 
*C>.  rinsed,  placed  la  a  physiological  salt 
da  to  30- mix  radiation.  AM  at  18  Ns.  in  8  TTM 
he  solution  wet  then  assayed  for  tedieactive 
entity  aeries  was  teeted.  Enhanced  afflux  of 
73  and  3.13  adi/an  ;  no  chants  wea  ebaarvad  at 
.•■xer- density  it  canverted  to  spacific  abaorp- 
*  ;  a  valid  mm h  too  bar  far  geoarsllrwi  haetu* 
•lying  coute  of  this  ratponaa.  Implication*  of 
adawa"  in  which  tha  offset  occurs,  both  tot  9/ 
IP*.  WMI,  •111  ho  diacuoaed.  Thla  la  an 
auilaa  and  data  not  necaxacrily  xaflact  EPA 


t-OI-A  mOBDUAVt'IXCUCXti  BESBQXTZ  2ATX0R  0J  AOETTLCaOLIItt  RKZFTgR 
CHANNEL  IP  CULTVR0  QUAIL  H70TUBC3.  P.  garnardi  ,  0.  D’lnaao  .  I.  Eueafci  . 
T  Areas!*.  C.  T^urslle  .  Department  of  Electronics ,  Uni  were  ity  of  E OM 
"la  Sapiens*",  Ml**  9mm. 

The  effects  ef  nicrowevea  an  the  acetylcholine  (ACh)  receptor  charatel  have 
been  atMiaed  veaerdina  ACh-lnduaed  currents  txperiesmts  war#  performed  on 
•  inale  3-3  day  aid  «a>stai*ted  ax>d  nen-imerveted  quail  nyotubes.  The 
properties  of  oinglo  ACh  activated  channel*  wore  investigated  using  the 
pstch-elssp  technique  with  high  resistance  oeale  betweai  the  naohria  and 
the  pipette  In  cell  attached  configuration,  aylgard  costed  pipettea  ware 
filled  with  scatylcholine  at  the  concentration  SxlO'8M.  Siaglo  ehatxiel  were 
analyied  by  the  threshold  crossing  Bathed.  Of  adenaaow  (7  Mf  at  10.73  CBr) 
ware  applied  using  a  bom  enterxia  (16  d*  gain)  at  a  distmee  of  23  wo.  from 
the  sample,  30*  iocllMted  to  the  ground  and  with  tha  C  vector  lying  on  the 
plane  of  Incidence.  In  the  control,  before  the  aorposure.  ACh-ectivatod 
thermal  currant  of  1.7*  +  0.01  pA  (ns so  ±  a.a.n. ;  7*0  evwits ) .  s  naan 
ehannal  open  tin*  of  3. *6  ±  0.13  an,  a  naan  closed  tine  of  18.7  +  1.)  ns 
wax*  neasured.  Mhen  nystuhea  were  exposed  to  Ml  for  120  o,  the  c Loved  tine 
increased  to  30. 1  ±  3.1  me  (300)  within  00s  fton  the  met.  Channel  open 
tine  (3.43  1  0  18  me)  and  current  (1.77  ±0.01  pA)  were  not  aKected  by 
irradiation.  After  the  and  of  exposure  a  partial  recovery  of  the  activity 
was  obtained.  It  la  concluded  that  wpoc*  to  adimnvoa  increases  the  mo 
of  deeaaultlaation  and  decreases  the  charnel  opening  probability.  Other 
asperlnantc  Indicate  that  thla  affect  develops  snip  within  range  of  field 
intensity 


l-aa-a  efpxcts  or  x  *3  ana  expobore  or  oportaeboub  ruman  immxiTl 
TRANSFORMATION  IN  VITRO.  E.C.  Elaon,  Walter  Kaed  AIR.  Silver  Spring.  IV 
20803.  E.M.  Ctcrakc,  C.C.  Devi*.  Elactrical  biglnacring  Dept..  University 
af  Maryland. .  Call***  Mrk.  MO  307*2.  P.  Csaraki.  CDM/FDA,  Rockville.  W 
20837 . 

Rntcn  peripheral  bleed  lymphocytes  in  ahranosonc  nedl«m  1A  without  MIA 
(Sibco)  hove  boon  exposed  to  continuous  (CM)  or  pulaod  (MO  2.43  CBs 
nicrowovas  using  the  autonatad  dealnctry  exposure  ay* tan  described  la  a 
companion  paper  (Joyner  at  *1..  thla  naetijg)  at  bestir*  md  nm-baaung  am 
levels  far  period*  up  to  120  hours.  At  harvest,  standardised  cytological 
preparations  wars  node  using  a  cytotpin  II  aytacantrifuge,  air-dried,  fixed 
in  noth  anal  and  stained  with  Wrlght/Olaaea  otaln.  The  percentages  ef 
witraMforsied  lymphocytes,  interned  lata  and  transfomed  Ipgmdblaetoid  eelLx 
wore  datamlned  baaed  on  narphologieal  features  quantitated  wslng  an 
"Opt  oner"  laaga  analyser  and  canparad  batwean  expoaed  <CW  or  |H)  end  than- 
exposed  tanplat.  CM  or  Ml  exposure  at  the  tana  average  SAR  resulting  In 
•  toiler  tanperatura  prof  ilea  over  tine  affect*  Ip^hocyt*  transformation 
differently.  Mechanistic  lapllcatiena  of  this  finding  are  discussed 
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tw-os  automatic  Doeitcntt  or  niapMve  ak>  imm  sxrosuu  or  sxotooi* 

CAL  SAMRXS  IN  VITkO.  K.N.  Joyner .  Tvltew.  Melbourne,  AiilnUa  B.C^ 
Elton.  Halt**  Boon  AIR.  Silver  Sprirw.  tO  IW05.  C.C.  p*vi».  S.H.  Cnnta, 
Eloctrleoi  Ent.  Dept..  »  of  ,  Cellaa*  Park,  K>  S07*2.  f.  ClIlUl, 
C£*fl/rDA.  Rockville.  IV  20857. 

A  nierovav*  avpeeur*  aid  dailoatlli  ajnta  will  he  dtsctlM  that  provide* 
non-porturbln*  autetttloO  recording  of  the  tRpmtm  hi » Lory  of  wpotad  end 
control  tablet  In  vitro.  Boeerda  of  tba  tmp*ratuj*«T>/tim(t)  behavior  of 
both  • ample  and  central  era  node  Wafer* .  Nurias,  aid  after  expeaur*.  f* cm 
the  recorded  T(t)  varaua  t  durin*  liaatla*  and  seeling  the  specific  absorp¬ 
tion  rat*  (SAB)  far  nltroweves  (or  tba  equivalent  energy  absorption  Nurias 
conventional  baatint)  ca  be  determined.  The  cepoaure  cystm  comet*  of  uaa 
aborted  aactiona  of  S  band  roctansular  wavesuld*  with  notched  coaxial  feed* 
placed  in  a  COj  tlseua  oulture  incubator  tharmoateted  at  }/.3'C.  Sabpla 
haldere  ere  eupported  centrally  in  the  wcvesulda  IB  blOCka  of  low-density 
expanded  polratyrena .  An  uacillator.  TWI  lapli fler .  isolator,  and  forward 
and  reflected  puwer  waters,  provide  the  power  to  tfto  epnart  wifiud*  Dio 
aaeand  wevoeuide  house*  tba  control  swnplaa .  Tanperaturos  m  ana  ample 
chamber  in  vevli  wava(ulda  ora  recorded  with  0  Nerda  Mdsl  HUB  iwrjetdt- 
in*  probe.  Tlio  tawparacure  XMdi«|i  are  dliltllM  and  atorad  in  a  desktop 
aenputar.  linear  least  aviate*  analysis  of  IttJ  warn*  t  with  the 

osemption  of  esnstant  power  absorption  end  Andorra  la>  of  mliif  pMM 
a  very  reliable  value  for  tho  SAR.  On-line  dotami nation  of  tho  SAB  in  this 
way  lUnn  accurate  assesmont  of  exposure  lailfonslty  Cron  one  ample 
chanhe?  to  the  ether,  and  permit*  accurate  calibration  of  tho  oyeten  In 
tame  of  e  fiAB/ferwasd-power  ratio. 
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Abstract — A  waveguide  exposure  system  with  automated  sample  temperature  measurement  it  described.  This 
system  provides  on-line  determination  of  the  temperature  profile  over  time  of  biological  samples  fa  Wire.  It  allows 
automated  computation  of  the  specific  absorption  rate  determined  from  beating /cooling  curves,  uses  mkaimally- 
perturbing  thermometry,  is  biocompatible  and  can  be  used  for  measurements  of  both  microwave  and  conventional 
beating. 


INTRODUCTION 

Dielectric  heating  is  a  recognized  mechanism  for  the 
induction  of  radiofrequency  (RF)  bioeffects  considered 
to  be  thermal  in  origin — that  is,  attributable  to  temper¬ 
ature  increases.  However,  apart  from  RF-heat  cell  killing 
(Sapareto  et  al.  1982;  Chang  et  al.  1987),  very  little  in¬ 
formation  is  available  on  quantitative  relationships  be¬ 
tween  RF  bioeffects  at  the  cellular  and  molecular  level 
and  the  temperature  profile  over  time  (TPT).  The  deter¬ 
mination  of  RF  heating/cooling  curves  can  be  used  to 
compute  the  specific  absorption  rate  (SAR  W  kg"')  and 
the  specific  absorption  (SA  J  kg"1)  in  an  exposed  sample 
(Stuchly  and  Stuchly  1986).  Biological  variables  and  ex¬ 
posure  conditions  can  be  controlled  in  experiments  with 
in  vitro  systems  to  a  degree  not  achievable  in  vivo.  An 
attractive  model  to  study  RF  bioeffects  is  the  transfor¬ 
mation  of  lymphocytes  in  vitro,  provided  the  biological 
variables  can  be  related  to  dosimetric  quantities  that 
characterize  exposure  (Czerski  1975;  Budd  and  Czerski 
1985 ).  To  accomplish  reliable  dosimetry,  we  designed  an 
exposure  system  with  provision  for  real-time  temperature 
monitoring  with  RF-field  non-perturbing  temperature 
probes.  The  exposure  system,  which  will  be  described  in 
detail  later,  has  multiple  sample  chambers.  To  allow  on¬ 
line  thermometry  and  dosimetry,  one  of  these  chambers 
is  used  as  a  site  for  a  non-perturbing  temperature  probe. 


From  the  temperature  ( T)/time(f )  history  of  the  sample 
chamber,  an  exposure  dosage  can  be  determined  and  de¬ 
scribed  in  terms  of  TPT,  SAR,  SA,  or  electric  field 
strength.  The  system  can  be  used  for  studies  of  RF  bioef¬ 
fects  in  any  tissue  culture  cell  line  or  other  in  vitro  bio¬ 
logical  samples.  The  temperature  monitoring  part  of  the 
system  and  the  associated  automated  data  analysis  can 
be  applied  to  the  study  of  temperature-dependent  effects, 
irrespective  of  the  modality  used  for  heating. 

EXPOSURE  SYSTEM 

A  schematic  diagram  of  the  system  for  microwave 
exposure  of  samples  is  shown  in  Fig.  1.  The  output  from 
a  CW  oscillator11  feeds  a  travelling-wave-tube  ( TWT )  am¬ 
plifier!  Pulsed  or  amplitude-modulated  waveforms  are 
obtained  through  the  use  of  a  p-i-n  modulator  and  driver#. 
For  high  peak  power,  but  low  duty-cycle,  pulsed  exposures 
the  signal  can  be  further  amplified  with  a  1  kW  ampli¬ 
fier*  •  The  amplified  signal  enters  a  shorted  section  of  S- 
band  rectangular  waveguide  300  mm  long  through  a 
matched  coaxial  feed.  The  waveguide  sample  holder  re¬ 
ceives  its  power  through  an  isolator  and  a  dual-directional 
coupler  that  allows  forward  and  reflected  powers  to  be 
measured  with  power  meterstt  end  waveforms  to  be 
monitored  with  a  diode  detectorft  and  oscilloscope.  A 
sample  holder  for  biological  specimens  is  supported  cen- 


*  Hewlett-Packard  Model  8616A,  Hewlett  Packard  Co..  3153  Porter 
Drive.  Palo  Alto.  CA  94304 

1  Hu|he*  Model  I I77H.  Hushes  Aircraft  Co.,  Electron  Dynamics 
Div.,  P.0  Boi  2999,  Torrance.  CA  90509. 


0  Hewlett-Packard  Model  540 JA. 

•*  MCL  Model  1 0214,  MCL  La  Granae.  IL  10284. 
ft  Hewlett-Packard  Model  432A. 

88  Hewlett-Packard  Model  423B. 
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Fig.  1.  Schematic  diagram  of  essentia]  components  of  microwave 
exposure  system. 


Microwave  Power  In 


Fig  2  Waveguide  exposure  assembly  showing  location  of  four- 
chamber  sample  holder. 


trally  in  the  waveguide  in  a  block  of  low-density  polysty¬ 
rene  that  fills  the  cross-section  of  the  guide,  as  shown  in 
Fig.  2.  An  identical  waveguide  section  with  an  identical 
plastic  sample  holder  is  used  for  sham  exposures.  During 
microwave  exposures,  both  waveguide  assemblies  were 


H  Mil«t  Scientific  a  4804.  Miles  Scientific.  Division  of  Miles  Lab¬ 
oratory  Inc..  Naperville.  IL  60566 

'  The  Narda  Microwave  Corporation.  435  Moreland  Road.  Haup- 
pa*e.  NY  11788 
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City.  UT  14108 


March  1989,  Volume  56.  Number  3 


Fig.  3.  Schematic  diagram  of  four-chamber  sample  holder  with 
minimally-perturbing  temperature  measurement  and  automated 
dosimetry  system. 


housed  in  a  C02  tissue  culture  incubator  thermostated  at 
37.0°C.  Various  sample  holders  were  tested:  a  four- 
chamber  plastic  tissue  culture  slide§§,  shown  in  detail  in 
Fig.  3,  was  found  satisfactory  in  our  application.  Each 
chamber  holds  1  mL.  Since  this  sample  holder  occupies 
only  about  10%  of  the  cross-sectional  area  of  the  wave¬ 
guide,  field  uniformity  from  one  chamber  of  the  sample 
holder  to  the  other  should  be  within  about  ±6%  of  the 
average  and  power  densities  should  be  within  ±12%  of 
the  average.  This  expected  exposure  uniformity  was  borne 
out  by  quantitative  dosimetry.  A  miniature  thermistor 
temperature  probe  enters  each  waveguide  section  through 
a  small  hole  in  the  top  shorting  plate.  The  hole  is  drilled 
nes&the  guide  wall  at  the  center  of  the  shorter  dimension. 
This  is  a  low  electric  field  point,  and  the  hole  produces 
minima]  disturbance.  The  active  end  of  the  thermistor 
probe  enters  the  sample  in  one  of  the  four  chambers  of 
the  tissue  culture  dish.  We  could  move  the  probe  from 
chamber  to  chamber  to  check  exposure  uniformity.  Two 
types  of  minimally-field-perturbing  probe  have  been  used 
in  this  way:  a  Narda  Model  801  IB  non-perturbing  double 
temperature  probe""  or  two  Vitek  Model  101  probestH. 
The  temperatures  in  both  the  exposed  and  sham-exposed 
samples  are  recorded  continuously  during  an  experimental 
run.  The  temperature  probes  are  connected  to  a  relay  ac¬ 
tivator##  which  is  connected  to  a  digital  voltmeter 
(DVM  )•  *  •.  Both  the  relay  activator  and  DVM  are  under 
the  control  of  a  desk-top  computerttt  through  the  IEEE- 
488  bus.  Under  computer  control,  the  temperatures  of 
exposed  and  sham-exposed  samples  are  recorded  sequen¬ 
tially  and  itored  in  memory.  Since  the  temperature  of 
sham-exposed  samples  is  essentially  constant,  it  is  mon- 


#*  Hewlett-Packard  Model  59306V 

•••  Kehhley  Model  192.  KeitMey  Instrument*.  Inc..  28775  Aurora 
Road.  Cleveland.  OH  44139. 

ttt  Hewlett-Packard  Model  16. 
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itored  less  frequently  than  the  temperature  of  the  exposed 
samples.  A  typical  experimental  protocol  would  involve 
10  temperature  readings  of  the  exposed  sample  for  every 
one  reading  of  the  sham-exposed.  For  conventional  ther¬ 
mal  exposures,  the  temperature  can  be  elevated  by  in¬ 
creasing  the  temperature  of  the  incubator. 

TEMPERATURE /TIME  ANALYSIS 

During  a  typical  experiment,  the  temperature  is  re¬ 
corded  at  regular  intervals,  with  a  minimum  measurement 
interval  of  <1  s  before,  during  and  after  the  exposure. 
The  “on”  and  “off”  times  of  exposure  are  recorded  on 
the  computer  by  the  operator  using  a  “soft-key"  interrupt 
capability.  At  the  conclusion  of  a  run,  the  7"(r)  behavior 
is  analyzed  to  determine  SAR. 

If  the  SAR  is  S  ( W  kg'1)  and  the  specific  heat  of  the 
sample  is  C  (J  kg'1  K'1),  then  the  rate  of  heating  during 
microwave  exposure,  for  small  temperature  variations 
over  which  C  can  be  assumed  to  be  constant,  is 


If,  in  the  absence  of  deliberate  heating,  a  sample  is  above 
or  below  its  equilibrium  temperature,  with  its  sur¬ 


roundings,  then  for  small  temperature  differences  from 
equilibrium,  its  natural  rate  of  temperature  change,  from 
Newton’s  law  of  cooling,  is  of  the  form 

(? m 

where  a  is  a  constant  that  depends  on  sample  geometry, 
insulation,  etc.  In  general,  when  microwave  power  is  ap¬ 
plied  to  a  sample  that  was  not  in  thermal  equilibrium 
with  its  surroundings  at  the  start  of  exposure: 

dT  S 

^mC  +  a(T«-T).  (3) 

The  solution  to  eqn  ( 3 )  can  be  written  in  the  form 

T-ro-^+r^-rojd -c-'j,  (4) 

where  T0  is  the  control  temperature  at  r  *  0.  Either  eqn 
(3)  or  eqn  (4)  can  be  used  to  determine  the  SAR.  From 
eon  (3)  we  have 


X 


O 


a 

E 

(S 


d  —  -i  <\i  »\i  ri 


Time  (s)  x  100 

Fig.  4.  Temperature/time  profile  of  a  sample  that  was  subjected  to  RF  heating  (solid  curve).  The  dashed  curves 
are  least-squares  linear  fit*  to  the  experimental  data  To  the  left  of  the  time  at  which  RF  power  was  applied,  the 
sample  was  cooling.  The  change  in  slope  at  the  initiation  of  the  RF  heating  provides  the  quantitative  dosimetry. 
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where  subscript  S  indicates  the  application  of  microwave 
power  corresponding  to  SAR  5.  Therefore,  determination 
of  the  rate  of  temperature  change  before  and  during  ap¬ 
plication  of  microwave  power,  or  during  and  after  the 
application  of  microwave  power,  determines  the  SAR. 

If  the  sample  is  in  thermal  equilibrium  before  ex¬ 
posure  begins,  eqn  (4)  reduces  to  the  simpler  form 

T-To-^d-e-'),  (6) 

which  can  be  used  to  find  S'  from  T(t).  In  general,  the 
approach  using  eqn  ( 5 )  is  preferable,  since  for  times  near 
/  *  0,  the  exponential  can  be  expanded  in  quadratic  form. 

In  practice,  this  analysis  is  performed  automatically. 
At  the  end  of  a  run,  the  point  of  inflection,  or  turning 
point,  in  T(t)  is  found  numerically.  This  is  done  by 
checking  the  change  in  slope  of  the  temperature  as  a  func¬ 
tion  of  time.  In  practice,  this  change  in  slope  occurs  im¬ 
mediately  upon  the  application  or  removal  of  microwave 
power,  within  the  time  resolution  of  the  data  acquisition 
system.  The  temperature  profile  to  the  left  and  right  of 
this  point  is  fitted  by  least-squares  to  a  linear  or  quadratic 
function,  and  the  change  in  slope  at  the  beginning  of  the 
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exposure  period  yields  the  SAR.  Figures  4  and  5  are  ex¬ 
amples  of  such  a  procedure.  Figure  4  shows  the  temper¬ 
ature/time  history  of  a  sample  that  was  cooling  prior  to 
the  start  of  microwave  exposure.  The  linear  fits  to  the 
cooling  and  h  ing  portions  of  the  curve  yield  the  SAR. 
Figure  5  she  *  ■  .he  heating  of  a  sample  that  was  quite  well 
equilibrated  oefore  exposure.  Because  the  signal/noise 
ratio  is  high,  the  SAR  was  determined  by  a  quadratic  fit 
to  the  heating  portion  of  the  curve.  The  quality  of  the  fit 
near  the  turning  point  can  always  be  examined  to  deter¬ 
mine  whether  the  slope  value  is  realistic.  However,  unless 
the  temperature /time  profile  has  very  low  noise,  it  is  gen¬ 
erally  better  to  use  the  linear  fit  to  7X0  near  the  turning 
point,  otherwise  a  biased  estimate  of  the  slope  can  result. 
Figure  5  also  shows  the  temperature/time  history  of  the 
sham-exposed  sample  to  illustrate  its  temperature  stability. 

RESULTS  AND  DISCUSSION 

The  advantages  of  the  procedure  described  above  are 
severalfold.  The  sample  need  not  be  equilibrated  before 
SAR  is  determined.  Repeated  determination  of  SAR  can 
easily  and  quickly  be  made,  which  allows  the  SAR  uni¬ 
formity  from  one  sample  chamber  to  another  be  deter- 
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Fig.  5.  Temperature /time  profile  of  expoied  (upper  solid  curve)  and  sham-exposed  samples  (lower  dashed 
curve).  The  dosimetry  for  the  exposed  sample  has  been  determined  from  a  quadratic  fit  to  the  heating  curve 
during  microwave  exposure  and  a  linear  fit  to  the  equilibrated  portion  of  the  curve  prior  to  exposure  (dashed 

upper  curve). 
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mined.  Exact  knowledge  of  the  microwave  power  is  not 
required;  reproducible  exposures  can  be  made  at  known 
SAR  by  using  the  measured  forward  power  corresponding 
to  a  given  SAR  measurement,  provided  the  experimental 
arrangement  is  not  altered  between  exposures.  Dosimetry 
is  not  affected  by  other  losses  in  the  system.  The  waveguide 
exposure  system  with  four  sample  chambers  was  found 
to  give  an  SAR  uniformity  from  one  chamber  to  the  other 
within  10%.  Thus,  experimental  samples  can  be  exposed 
and  examined  in  triplicate,  the  fourth  chamber  being  used 
for  the  insertion  of  the  temperature  probe. 

The  system  is  biocompatible  and,  depending  on  cell 
line,  cell  density  and  medium,  allows  continuous  expo¬ 
sures  of  several  days  duration.  The  biocompatibility  of 
this  arrangement  was  proven  by  studying  the  growth  of 
human  lymphocytes  under  various  conditions  within  the 
waveguide.  A  concentration  of  106  cell  mL"'  in  chro¬ 
mosome  medium  lA$t$  with  or  without  phytohemag- 
glutin  was  incubated  for  72  or  120  h  in  the  exposure  sys¬ 
tem  with  no  applied  RF  power.  Cell  morphology,  lym- 
phoblastoid  transformation,  mitotic  index  and  cell 
viability  tested  by  neutral  red,  Janus  green  and  trypan 
blue  stains  did  not  differ  from  those  in  lymphocyte  cul¬ 
tures  from  the  same  donor  cultured  in  a  conventional 
CO;  incubator  according  to  the  protocol  provided  by 
Gibco  with  their  media.  As  well  as  providing  biocompat¬ 
ible  exposure  conditions,  the  system  is  flexible.  It  allows 
the  exposure  of  tissue  culture  cells  growing  in  suspension 
or  in  monolayers  and  can  easily  be  used  with  cell-free 
samples. 

The  exposure  system  is  relatively  simple  and  can  be 
assembled  from  off-the-shelf  components.  Its  capabilities 
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can  be  expanded  by  introducing  computer  control  of  RF 
power  input  based  on  feedback  from  temperature  mea¬ 
surements  in  a  similar  way  to  the  system  described  by 
Chang  et  al  ( 1 987 ) .  This  system  used  a  stripline  operating 
at  915  MHz  and  allowed  uniform  exposure  of  10  mL 
tissue  culture  samples.  The  temperature /time  profile  of 
heating  could  be  pre-programmed,  and  biocompatibility 
allowed  exposures  of  several  days  duration.  A.  W.  Guy 
(1977)  designed  a  system  for  exposure  of  5  mL  samples 
in  the  range  from  dc  to  I  GHz.  Field  strength  dosimetry 
was  obtained  from  on-line  measurement  of  feed-line 
impedance.  The  temperature  of  the  sample  could  be 
maintained  at  a  preset  level  using  a  circulating  liquid  heat 
exchanger.  No  data  on  biocompatibility  of  this  system 
were  provided,  however  the  information  given  suggests 
suitability  only  for  short-term  exposures.  J.  C.  Lin  (1976) 
and  Chen  and  Lin  (1978)  proposed  the  use  of  micropi- 
pettes  immersed  in  a  fluid-filled  waveguide  irradiation 
chamber  for  in  vitro  studies  under  controlled  temperature 
and  dosimetric  conditions.  This  system  was  also  only 
suitable  for  short-term  exposures.  J.  W.  Allis  et  al.  (1975 ) 
installed  a  waveguide  exposure  system  in  a  Cary  15  spec¬ 
trophotometer  allowing  simultaneous  exposure  and  spec¬ 
trophotometric  biochemical  determinations.  Our  system, 
as  well  as  those  described  in  the  papers  quoted  here,  pro¬ 
vides  a  wide  range  of  possibilities  for  in  vitro  studies  under 
controlled  conditions,  including  precise  RF  dosimetry. 
By  using  one  of  these  systems,  or  combining  desired  fea¬ 
tures  of  some  of  them,  most  requirements  for  the  study 
of  RF  bioeffects  at  the  cellular  and  /  or  molecular  level  can 
be  met. 
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